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1.  INTRODUCTION 

Lattice  materials  are  porous  materials  consisting  of  periodic  cells  or  non¬ 
periodic  cells.  The  cells  are  composed  of  rods,  or  shells,  or  solid  structures. 
The  size  of  cell  is  usually  small  with  respect  to  the  size  of  the  body  filled 
with  the  lattice  materials.  The  lattice  materials  with  simple  microstructures 
are  characterized  by  a  single  length  scale,  for  instance,  Lattice  Block  Ma¬ 
terials  which  are  developed  by  JAMCORP  corporation.  The  hierarchic  lat¬ 
tice  materials  have  hierarchic  multi-scale  microstructures.  In  either  case,  we 
deal  with  a  multi-scale  problem.  The  lattice  materials  can  offer  significantly 
higher  strength-to-weight  and  stiffness-to-weight  ratio  than  their  base  materi¬ 
als  and  have  obvious  advantages  in  engineering  applications.  Besides  heteroge¬ 
neous  materials[17,  28,  30,  33],  lattice  models  are  used  in  many  other  applica¬ 
tions  such  as  porous  media[12,  16,  33],  fracture  models[31],  crystal  physics[21], 
biophysics[20].  For  a  survey  of  some  applications,  we  refer  to  [30,  32]. 

Various  mechanical  models  for  the  lattice  problems  with  periodic  and  non¬ 
periodic  structures  have  been  studied.  There  are  papers  addressing  these  prob¬ 
lems,  especially  in  the  mechanics,  material  science,  and  physics  literatures.  For 
mathematical  theory  which  is  related  to  the  problem  of  the  lattice  materials 
we  refer  to  the  book  [8]  and  her  various  papers,  e.g.  [7,  9,  10].  Various  math¬ 
ematical  aspects  of  lattice  problems  have  been  studied.  Problems  of  periodic 
lattices  in  entire  spaces  and  numerical  methods  are  analyzed  based  on  Greens 
function  and  Fourier  transform  [3,  21,  22,  23,  24].  A  combination  of  homoge¬ 
nization  and  multigrid  method  on  unstructured  mesh  is  applied  to  problems  of 
non-periodic  lattices  [6, 14].  The  recovery  method  has  been  developed  recently 
for  problems  of  bounded  lattice  [4],  which  convert  discrete  lattice  equations  to 
a  finite  element  discretization  of  a  continuous  partial  differential  equation. 

The  results  of  these  papers  are  mostly  for  problems  of  periodic  lattices  in 
entire  spaces  Rd,  but  not  for  problems  of  unstructured  lattices  on  unbounded 
and  bounded  domains  with  prescribed  boundary  conditions.  In  addition,  the 
scale  of  cells  in  these  papers  is  assumed  so  small  that  asymptotic  arguments 
such  as  homogenization  can  be  utilized.  Further,  most  papers  only  address 
the  problems  in  presence  of  absolute  terms  in  the  equations  so  that  the  cor¬ 
responding  bilinear  form  satisfies  the  inf-sup  condition  on  a  pair  of  discrete 
Sobolev  spaces.  In  practical  applications,  these  assumptions  may  not  be  valid. 
If  the  scale  of  cells  may  not  extremely  small,  then  the  homogenization  will  not 
give  us  satisfactory  solution  of  lattice  problems.  In  particular,  the  absolute 
terms  do  not  exist  in  applications,  it  is  a  pure  mathematical  hypothesis  for 
the  convenience  in  mathematical  analysis.  For  problems  without  such  abso¬ 
lute  terms  we  need  to  develop  new  mathematical  approaches  for  problem  of 
unstructured  lattices. 

This  final  progress  report  summarizes  the  progresses  of  our  research  in  the 
period  1999-2002,  including  many  results  which  are  original  and  have  not  been 
published  yet.  In  the  past  three  years,  we  started  with  problems  of  periodic 
lattices  in  entire  spaces,  and  then  focused  on  the  mathematical  framework 
for  problems  of  unstructured  lattice  in  unbounded  domains  without  absolute 
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terms.  The  Fourier  transform  is  a  very  powerful  tool  for  study  of  problems  of 
lattices  with  periodic  structure  and  for  design  of  effective  numerical  methods 
such  as  generalized  p- version  of  finite  element  methods  [25,  26] ,  even  for  those 
on  bounded  and  unbounded  domains  with  boundary  layers.  But  it  is  not  valid 
for  unstructured  lattices.  The  new  approach  is  to  establish  the  equivalence 
between  the  problems  of  lattices  and  partial  differential  equations  (not  homog¬ 
enized  equations),  and  the  equivalence  between  grid  functions  and  continuous 
functions.  With  help  of  such  equivalences  we  can  establish  various  embedding 
results  among  discrete  function  spaces  which  lead  to  the  existence  and  unique¬ 
ness  of  solution  of  lattice  problems,  also  it  will  results  in  an  effective  algorithm 
via  the  numerical  solution  of  the  corresponding  partial  differential  equations, 
e.g.  multigrid  method.  The  problems  we  presented  are  of  truss  type  which  re¬ 
sult  in  systems  of  difference  equations  with  infinite  number  of  unknowns.  Such 
lattice  problems  are  addressed  in  discrete  energy  space  or  Sobolev  spaces  over 
lattices.  The  mathematical  framework  we  have  established  and  mathematical 
foundation  we  have  laid  down  for  the  problems  of  unstructured  trusses  can  be 
applied  to  or  generalized  to  problems  with  complicated  microstructures  such  as 
plates  and  shells,  or  three  dimensional  solid  structures  and  many  non-periodic 
lattice  problems. 

The  report  is  organized  as  the  follows.  In  Section  2,  we  address  problems 
of  periodic  trusses  in  entire  spaces  Rd,  d  —  1, 2, 3.  The  problem  is  setup  with 
various  important  concepts  such  as  connectivity  and  rigidity.  With  Fourier 
transformation  the  equivalence  between  a  discrete  problem  and  a  semi-discrete 
problem  is  established,  and  a  representation  formula  of  solution  of  truss  prob¬ 
lem  is  derived,  at  the  end  of  this  section,  two  lattice  problems  are  analyzed 
in  the  theoretical  framework.  In  next  section  we  present  the  major  progresses 
of  research  on  the  problems  of  unstructured  trusses  without  absolute  terms  in 
entire  spaces.  The  extension  of  grid  functions  to  continuous  functions  over  Rd 
is  defined  by  a  linear  interpolation  is  explicitly  constructed,  which  leads  to  an 
embedding  of  the  energy  space  into  weighted  L 2  spaces.  The  existence  and 
uniqueness  of  the  solution  are  proved  for  problem  without  absolute  terms  in 
one,  two  and  three  dimensions  by  different  way.  In  Section  4  a  boundary  value 
problem  of  unstructured  trusses  in  half  spaces  R+  is  analyzed.  A  sufficient 
condition  on  the  external  force  /  is  derived,  under  which  the  existence  and 
uniqueness  of  solution  of  the  boundary  value  problem  is  proved.  At  last  we 
conclude  with  current  and  future  research  directions  on  problems  of  lattices, 
mathematical  models  and  computational  methods.  In  two  appendices,  the  so¬ 
lution  of  Possion  equation  in  entire  spaces  Rd,  d  =  2, 3  and  the  boundary  value 
problem  in  half  spaces  Rrf,  d  —  2,3  are  discussed  in  modified  function  spaces, 
which  are  parallel  to  the  problems  of  lattices  in  entire  and  half  spaces,  some 
embedding  results  are  essential  to  the  existence  of  the  solution  of  such  lattice 
problems. 
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2.  PERIODIC  LATTICES  IN  ENTIRE  SPACES 


2.1.  General  setting  of  periodic  in  entire  spaces  Rrf 

Lattices  are  comprised  of  cells  and  nodes.  Let  Q  be  a  master  cell  in  Rd 
with  unit  size,  which  is  an  interval  in  one  dimension,  and  is  a  polygon  in 
two  dimensions,  and  is  an  polyhedron  in  three  dimensions.  The  master  cell  is 
extended  periodically  in  entire  space  by  an  integer  translation: 

d 

Qm  =  {y€Rd\y  =  x  +  ^2mit{l\xeQ},meZd.  (2.1) 

i—1 

where  Z  =  {0,  ±1,±2, . . . }  denotes  the  set  of  all  integers,  and  fW  is  a  unit 
vector  in  Xi  axis.  There  is  a  set  Kq  of  nodes  {x^}qK=l  in  the  master  cell  Q, 
and  a  set  Km  is  the  integer  translation  of  Kq  by 

d 

Km  =  {x^  =  xiK)  +  €KQ},m€Zd.  (2.2) 

i—  1 

Note  that  the  indices  k  of  nodes  in  each  cell  Qm  are  the  same  although  the 
locations  of  these  points  in  different  cells  are  different.  Hence  we  denote  the 
set  of  indices  {1,2, ...  ,q}  by  /C.  Without  losing  generality,  we  assume  that 
the  cells  Qm’s  and  sets  Km's  are  mutually  disjoint,  namely, 

Qn  C  Qm  =  0?  Kn  PI  Km  —  0,  Tl  ^  771,  77,  777  £  Z  . 

The  lattices  shown  in  Fig.2.1  and  2.2  are  two  typical  examples  with  Q  = 
[0, 1),  fC  =  {1, 2}  and  Q  =  [0,  l)2, 1C  =  {1},  respectively. 


E, 


E2 


(a)  A  lattice 


-W,  ■  f 

i-mP  f 


(b)  Master  cell  Q  furnished  with  springs 

Fig.  2.1  A  periodic  lattice  in  R1 


We  now  further  specify  the  connectivity  of  lattices.  By  b(-m'K,n,)'\we  denote 
an  elastic  rod  connecting  the  nodes  x^m,K^  and  x^n,x^  with  intersect  area  A  and 
length  By  E(-K,n,x^  we  denote  the  Young’s  modules  of  the  elastic  rod 

b(K,n,A).  por  the  8^  0f  simplicity,  we  will  omit  cell  index  m  whenever  777  =  0. 
For  instance,  we  write  =  x^  ,  b^0,K,n’^  =  b^K,n’^,  etc.  We  assume  that 

(C.l)  Each  node  is  connected  to  others  by  the  rods, at  least  one  node  and  at 
most  M  nodes; 

(C.2)  Any  two  nodes  x<'m,K'>  and  x('n,x')  are  linked  by  the  shortest  chain  Lm<K>n> 

x(m,K)  _  x(ni,\\)  x(n2M)  x(naM)  .  . .  — >  /p(»«iA»)  =  x(n>P  such  that 
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x(nt,\t)  is  connected  to  a;(n*+1>A*+1))  1  <  t  <  s  —  1,  and 

|  j /j»(^t+i >^t+i)  X(^)|, 

l<t<s-l 

where  rj  is  independent  of  m,  n,  k  and  A; 

(C.3)  The  length  of  rods  are  uniformly  bounded  i.e.  for  any  and  x^n,x^ 
which  are  connected,  there  holds 

bx  <  b^m’K’n’^  =  \xM  -  z(n’A)|  <  b2. 

To  effectively  describe  the  connectivity  of  the  lattice,  we  introduce  BK  and 
Bk,  a 

Bk  =  {(n,  A)  G  Zd  x  K  such  that  and  x^n,x^  are  connected  }  (2.3a) 

and 

Bk,\  =  {n  €  Zd  such  that  (n,  A)  G  BK}.  (2.3 b) 


(a)  A  lattice 


(b)  Master  cell  Q  furnished  with  springs 


Fig.  2.2  A  periodic  lattice  in  R2 


Bk  and  BKi\,  based  on  the  connectivity  of  the  nodes  x^0,K^  in  the  cell  Qo,  can 
be  periodically  generalized  to  sets  B^  and  B^  for  all  m  G  Zd  by  the  integer 
translation.  Due  to  the  periodicity  ,  it  is  easy  to  verify  that 

n  G  BKtx  if  and  only  if  -  n  G  BX,K  (2.4) 

and 

E(K,n  A)  _  £l(A,-7l,K)_  (2.5) 

We  further  assume  that  lattices  are  rigid  in  the  sense  that  if  the  new  config¬ 
urations  of  all  cells  resulted  from  a  continuous  transformation,  which  remains 
the  lengthes  of  all  rods  fixed,  are  congruent  to  the  original  ones.  A  cell  is  rigid 
if  the  new  configuration  of  the  cell  resulted  from  a  continuous  transformation, 
which  remains  the  lengthes  of  all  rods  of  the  cell  unchanged,  are  congruent  to 
the  original  one.  These  difinitons  of  rigidity  are  based  on  the  graph  theory  and 
coincide  with  those  of  [1].  Obviously,  a  lattice  is  rgid  if  all  cells  of  the  lattice 
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are  rigid,  but  a  rigid  lattice  may  have  some  non-rigid  cells.  For  instance,  a  cell 
consisting  of  triangles  in  two  dimensions  and  tetrahedrons  in  three  dimensions 
for  which  each  edge  is  an  elastic  rod,  is  rigid,  a  lattice  which  is  composed  of 
such  cells  is  referred  as  a  triangular  or  tetrahedral  lattice.  Of  course,  a  rigid 
lattice  may  not  be  triangular  or  tetrahedral.  For  the  various  definition  of  rigid¬ 
ity  in  general  and  verification  of  rigidity  for  a  lattice  or  a  truss  system,  which 
is  not  trivial,  we  refer  to  [2]  and  also  to  [1,  11,  27,  30]. 

A  lattice  is  characterize  by  the  local  structure  JC,  the  global  and  periodical 
translation  on  Zd,  and  the  connectivity  BK.  We  now  denote  the  lattice  with 
the  above  structures  by  G  =  G(JC,  2d,  BK). 

2.2.  A  Truss  Problem  on  Periodic  Lattice 

Let  u  =  (um)meZd  and  um  =  (um>K)Ke)c  be  a  grid  functions  on  Q  and  fC,  re¬ 
spectively,  and  each  um^K  is  a  vector  (ulm  ^  . . . ,  usm  K)T  with  s-components. 

In  one  dimension,  s  =  d  =  1,  and  umjK  denotes  the  displacement  for  elastic 
rods  or  the  temperature  of  heat  problems  at  the  node  x^m,K^ .  In  two  and  three 
dimensions,  s  =  1  when  um>K  denotes  the  temperature.  For  the  elastic  prob¬ 
lem,  s  >  d.  If  the  connections  of  rods  are  non-rigid,  then  s  =  d,  and  umA 
denotes  the  displacement  at  the  node  If  the  connection  of  rods  are 

rigid,  then  s  =  d(d+  l)/2  for  d  =  2,3,  and  um^  denotes  the  displacement 
and  rotation  at  the  node  x^m'K\  We  furnish  the  rods  with  springs  in  the  axis 
directions  at  each  node  with  Hook’s  coefficients  denoted  by  diagonal  matrices 
C<(m,/c)  _  q(k),  m  e  zd,  k  e  JC.  We  assume  that  the  ratio  of  the  length  of  rods 
and  the  intersect  area  A  of  rods  >>  1.  For  the  convenience  to  characterize  the 
nature  of  our  methodology,  we  will  focus  on  the  case  that  s=d  ,  namely,  the 
rods  join  with  hinges,  and  the  bending  is  not  considered  here. 


If  external  forces  exert  on  the  rods  at  the  nodes,  denoted  by  /  =  ( fm )  = 
(/m,«)(m,K)e^xX,  we  have  the  equilibrium  equation 

-  £  +°W^  =  /^.Vm  €  Z*,Vk  e  K. 


(n,X)€BK 


(2.6) 


with 


=  ,!£.(«...»  ~  *W)  ~  *W>T  (2.7a) 

|3j(n,A)  __  |a;(n’A)  — 

which  is  a  matrix  for  s  =  d  >  1  and  a  scalar  quantity  AE^K,n,>^  for  s  =  1, 1  < 
d  <  3,  and 

C(k)  =  diag(C[K),  Cf}, . . . ,  cj°),  Cf}  >  0.  (2.7 b) 

Similarly,  is  a  matrix  for  s  —  d  >  1  and  a  scalar  quantity  for 
s  =  1,1  <  d  <  3. 

Let  Hl(Q)  and  L2{Q)  be  the  Sobolev  spaces  over  the  lattice  Q  with  the 
norms 

imiLc)  =  £  £  K«i2  (2-8) 

m€Zd  «€/C 


8 


and 

IMIjr^c)  =  Mlr1^)  +  llulll2(e) 

where  \u\h^(q)  is  the  semi-norm, 

I  I2  —  V'  V'  l^m+n,A  ~ 

\U\Hl{g)  ~  2-~!  Z-/  2-^  \x(m+n,\)  _  x(m,K,)\2' 
m£Zd  k€)C  (n,A)€J9« 


(2.9a) 


(2.96) 


The  corresponding  variational  problem  is  defined  as 

B{u,  v)  =  F(v)  (2.10) 

with  the  bilinear  form  on  Hl{Q)  x  Hl{Q) 

D/,.  _  V'  V  /V1  ^  /-plK.n.X)  (um+n,x  ~  um,n) 

Jj\U,V)  Z^/m€Zd  2-JKEK\2-/(n,X)£BK  2'  |^(m+n,A)  _ ^.(m5K)  |  ’ 

(2.11a) 

and  the  linear  functional  on  H 1  (Q) 

n»)  =  D  o>  e.nfo 

m€Zd  k€K 

where  ( x ,  y)  =  Y2j=i  xjVi  is  the  inner  product  of  two  vectors  in  Rrf,  and  |x|2  = 
(x,x). 

The  energy  of  the  elastic  rods  is 
G(u)  =  B(u,  u)  =  £m€2. 


x(m,K.)  | : 


(Um+n,\  Um ,k)  x  /£(«),,  u  U 

|a;(m+n,A)  _  x(m>K)  |  /  +  Um’K ’  S 1 


The  energy  space  denoted  by  E(Q )  is  the  family  of  all  grid  functions  u  on  Q 
with  finite  energy  G(u),  and  ||t/|j£;(^)=  G(«)1/2  is  referred  as  the  energy  norm 
of  u. 


Proposition  2.1  If  a  lattice  Q  is  rigid,  then 


ED  E  <E<“ 


,n,A) 


m£Zd  K&JC  (n,A)€B/c 


(^m+n,A  (^m+n, A  v  q 

|/^(m+n,A)  _  j^(m+n,A)  ^(m,/c)|' 


(2.12) 


if  and  only  if  a  is  a  rigid  body  motion. 
Proof  The  proof  is  technical,  we  refer  to  [2] . 


Triangular  and  tetrahedral  lattices  are  rigid,  but  a  rigd  lattice  may  not  be 
triangular  or  tetrahedral.  It  can  be  proved  that  the  lattice  shown  in  Fig.  2.3 
is  rigid,  but  is  not  of  triangular  type.  The  lattice  shown  in  Fig.  2.4  is  neither 
rigid  and  nor  of  triangular  type. 


/]\  /[\  /|\  /[\  /|\ 


Fig.  2.3  A  rigid  and  non-triangular  lattice  in  R2 


i _ i 


T  T  T  T  T  T  T  T 

(a)  A  non-rigid  lattice 

1  1  i  1  1  1  1  J 


(b)  A  non-rigid  body  motion  on  a  non-rigid  lattice 

Fig.  2.4  A  non-rigid  lattice  in  R2 


Lemma  2.2  The  bilinear  form  B  on  Hl(Q)  x  Hl{G)  given  in  (2.10)  is  contin 
uous,  and  it  is  coercive  if  >  0  for  all  k  €  K. 


Proof.  Note  that 
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where  Qe  is  the  angle  between  the  rod  b^re,n,A-)  and  the  x^-axis.  We  have  imme¬ 
diately 

\B(U,V)\  <  C\\u\\Hi(Q)\\v\\Hi(g) 

with 

C  =  Max  {  max  A  E('li,n,x\  max  C^}. 

n£BK}X,K,\£K  1  <t<d,n£K 

Hence,  the  bilinear  form  B  is  continuous.  Since  >  0  for  all  k  E  1C,  we 
have 

-  IMIW 

Here  and  hereafter  ”=”  means  equivalent  with  constants  independent  of  major 
subjects,  e.g.  the  functions  u.  For  s  =  1,  there  holds 


/pfK.n.Al  (um+n,\  c)  (^ro+n, A  %a)  \ 

'  |^.(m+n,A)  _  2.(m,K)|’  |^;(m+n,A)  _  j.(m,K)y 


*)|< 


A  i?(K,n,\) 

=  |I(„+„,»-I(m,.)|4l<x(’"+"A)  - 

A£(K,n,\) 

=  ^(m+n.A)  _  x(m,K)\2  \Um+n,X  ~  wm,« |2  >  di\u\2Hl(g) 

with  d\  =  minneBK  XtK,,\eic  AE(-K,n,x\  For  s  =  d  >  1, 

/£'(,<c.n.A~)  (Um+n,X  ~  ^m,n)  (Mm+ii, A  ~  Mm,n)  ' 

'  |^.(m+n,A)  _ |  ’  |^(m+n,A)  _ 

4  pi(K,n,\) 

=  |;r(ro+..»_I(m.»)|4K:C<m+"'")  -  “">+”■*  -  “”*■«)  I' 

=  |a.(m+n,A)  _  jpCm.Kjp  lMm+n’A  “  ^m.KpCOS2^^  >  dlMffi(0) 

where  <4  =  min<t6Jc,(»,A)eB1«  AE^’n’x)cos2(f)Ktnix.  By  the  rigidity  assumption,  u 
will  be  a  rotation  if  di  =  0,  which  is  not  in  L2(Q).  Hence,  d\  >  0  for  u  E  Hl{Q), 
and 

B(u,u )  >  di\u\2Hi(g)  +  >  d||«||/fi(0) 

with  d2  =  mini <t<d,K£K >  0.  and  d  =  min{dl,d2}.  □ 


(2.13) 


Theorem  2.3  Suppose  that  C M  ^  0  for  all  k,  E  1C.  Then,  for  any  /  E 
,  the  variational  equation  (2.7)  has  a  unique  solution  u  E  Hl(Q), 

and 

In  particular,  if  f  E  L2(Q),  there  holds 

IMIjtHS)  <  C'll/lli2(e)-  (2-14) 


Proof.  The  theorem  follows  from  the  previous  lemmas  and  Lax-Milgram 
Theorem.  □ 

Remark  2.1  The  condition  that  >  0  for  all  k  E  K  can  be  weakened  to 
C Y*  ^  0,  the  energy  norm  ||w||f!(g)  is  equivalent  to  the  norm  ||u||i/i(a)  due  to  the 

connectivity  assumption  [2].  If  CjK^  =  0,  Theorem  2.2  can  not  stand  because 
the  energy  space  E(Q)  is  not  equivalent  to  Hl{Q)  and  it  is  not  embedded  in 
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L2(G).  Because  of  the  importance  of  lattice  problems  without  absolute  terms 
(i.e.  C jK^  =  0)  in  applications  we  shall  elaborate  it  in  Section  3  and  4. 

The  relation  between  the  solution  of  the  equilibrium  equation  (2.6)  and  the 
solution  of  the  variational  equation  (2.10)  is  given  in  next  Theorem. 

Theorem  2.4  If  u  G  Hl(G)  is  the  solution  of  the  variational  equation  (2.10), 
then  it  satisfies  the  equilibrium  equation  (2.6).  Vice  versa,  If  u  G  H1(Q)  solves 
the  equilibrium  equation  (2.6),  it  satisfies  the  variational  equation  (2.10). 

Proof.  We  first  prove  that  the  solution  u  of  the  equilibrium  equation  (2.6) 
satisfies  the  variation  equation  (2.10).  For  v  G  Hl(G),  multiplying  (2.6)  with 
vmik  and  summarizing  with  respect  to  m  and  k,  we  have 


y-'  /p(/c,n,A)  Mm+n, A  ~  _ vm,K _ \ 

ZL/  '  J  '  |/p(m+n,A)  _  I’ 

m€Zd  k,A€/C  neB^x  1  11 


m£Zd  kE/C 


m€Zd  «€/C 


The  above  sums  exist  since  u,  v  G  Hl{G)-  Letting  m  +  n  =  m'  and  n  =  —  n', 
we  get 


E\  '  V  s  /-p(K.n.A)  ^m+n,A  ^m,K 
m€Zd  J  '  |/^(m+n,A)  _  /p(m,K)  |  ’ 

K,\€Kn£BKjX  1 


/p(m,K)| 5  |^(m+n,A)  _  ^(m,/c)| 


) 


y''  /p(/t,-ra.A)  um',A  'Mm'+n',*: 

m'£Zd  /  v  /  '  |T.(m',A)  _ +re',«)  I  ’ 


^m'+n'  ,i c 


'4-n.'  kM  ) 


k,A€£  -n'eBK, a 


|j.(m',A)  _ ^.(m'+n',^)  |  ’  |,j.(m',A)  _ ^.(m'+n',K)  | 


Due  to  the  properties  (2.4)  and  (2.5)  : 


g(re,-n',A)  _  g(A,n',(c) 

-n'  G  BKjX  iff  n'  6  £a,* 


there  hold 


V  /xp(K-n;,A)  um',\  umf+n’,K 

mlezd  Z-/«,AE/C  l^-n'eBKtX\  \x(m,iX)  __  1 5 


|^(m',A)  __  | 5  ^(m^A)  __  ^(m'+n'jK)  | 

V  /F(Ay>«)  Um’iX  Z  _ \ 

m(EZd  AE/C  Z^n'EBx,*'  |^(m',A)  __  ^(m'+n'j/c)  |  ’  |,£(m',A)  _  ^m'+ra',*)  | ' 

(2.15) 

Which  leads  to  the  (2.10)  immediately. 

We  now  show  that  the  variational  solution  u  G  Hl(Q)  solves  the  equilibrium 
equation  (2.10).  Let  v  G  Hl{Q)  be  such  that  vm>K  =  0  for  all  k  G  /C  except 


) 


12 


) 


k  —  1.  Then  the  variational  equation  leads  to 

Ey  / 1  pj(/t,n,l)  ^m+n,!  ~  Mm,/c  _ 

2-j  2—j  '2  U(m+n,l)  _  klm+n,1)  —  x(TO>K)  I 

mezd  «€/C  n€B*,i  1  11 

Er1  / 1  E(l.n.l)  Mm+n, A  ~  ^m,l _ ^Vl _ \ 

/  ^  ^  '2  I I  ’  —  ^0^,1)  1 

mezd\eJCneB1<x  1  11  1 

d"  ^  '  (C^  ,  Urn,l)  =  ^  ]  (/m,l)  ^m,l) 

m€Zd  m€Zd 

Selecting  um>K  such  that  uTO)K  =  0  for  all  k  G  /C,  m  E  2  except  %,i,  Due  to 
(2.15),  we  obtain 

EY^  /■pOc.n.A')  Mm+w,A  ~  «m,l _ ^L1 _ \ 

/  j  '  |a>(m+n,A)  _  x(m,l)  |  ’  |a;(7n+n,A)  _ £;(m,l)  | ' 

Tl€:Brh,\ 

+  (C^  )  ^m,l)  ~ 

which  implies 

EV  |?(tt,n,A)  nm+n,A  ~  Mm,l _ l-C^U-  ,  =  f-  , 

U(m+n,A)  _  12  m>1  •”71’1 

ASB«  neBm.A  1  1 

Similarly,  there  holds  for  any  m  G  Zd 


_  £  £  +  CO^,,  = 


A€Bk  n£B1A 


Actually,  the  above  argument  can  be  carried  for  any  k  €  JC.  Thus,  we  have 
the  equation  (2.6).  □ 


2.3  Fourier  transform  for  lattice  problems 

For  grid  functions  u  on  the  periodic  lattice  Q  in  entire  spaces,  we  introduce 
the  Fourier  transform 

T{u)  =  J2  u™ei{m,t)  =  Ht)  (2.16a) 

m€Zd 

which  is  a  linear  functional  over  the  space  C™r(Id)  =  {u  G  C°°(Id)\u(t)  is  a 
27r-periodic  function  },  where  Id  =  (— n, 7r)rf.  u(t)  is  a  complex-valued  vector 
function,  (uK(t))Kefc,  and  each  uK(t)  has  s  components  ueK(t),  1  <  t  <  s,  and 


fi ,(t)  =  Y,  Vk  S  AC. 

mezd 

(2.166) 

The  inverse  Fourier  transform  is  defined  as 

=  ( um)meZd 

(2.17a) 

for  any  u{t)  G  C^r(Id),  and 

Um  =  (27r)"d  [  uity-^'^dt 

Jid 

(2.176) 
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For  lattice  problems  we  are  interested  in  some  of  specific  spaces  over  Q, 
e.g.  L2(G),  then  F(u)  6  L2(Id ),  which  has  a  stronger  topology  than  the  linear 
functional  on  the  space  C™r(Id).  In  particular,  we  are  interested  in  the  Fourier 
transform  on  the  spaces  L2(Q),v>  0  with  the  norm 

IMIl2(0)  =  +  \m\2Y\um,K\2. 

mezd  *e>c 


Lemma  2.5  The  Fourier  transform  realizes  isomorphism:  L2(Q)  -B-  L2(Id )  and 
L2U(Q)  b-  Hper(Id),  where  Hper(Id)  is  the  subspace  of  2ir-  periodic  functions  in 
Hv{Id),  and 

li2(S)  =  (^7r)  dll*llL2(/d)’ 


u 


n,  2  rsj  -  2 

U  LUG)  —  U  H^{Idy 


(2.18a) 

(2.186) 


Proof.  It  is  easy  to  verify  that 

iniV)  =  [  EEE Ki2*  =  /„£££ i 

v  '  J  Id  mezd  K,eic  o<e<s  JI  mezd  *e>co<e<s 

=  (27r)d  £  £  £  K,P  = 

m£Zd  kG/C  0<^<s 

which  implies  an  isomorphism  :  L2(£)  O  L2(Id)  and  (2.18a).  For  u(t)  6 
Hper(Id )  with  integer  v  >  0,  there  holds 

d 

Dau(t)  —  PJ(mi£)a*u(t) 

l-i 

for  any  a  =  (a-i,  a2,...,  ctd)  with  |a|  =  Y^i<i<d  a*  -  u->  which  leads  to 


i^ii  H^id)  —  ^2  + \m\2y\ 

mezd 


'U'm,  — 


iw 


For  non-integer  v  >  0,  H%er(Id )  is  defined  as  an  interpolation  space,  and  (2.18b) 
stands  for  non-integer  v  as  well.  □ 

We  now  apply  Fourier  transform  to  the  variational  problem  (2.10).  We 
introduce  a  bilinear  form  B  and  a  linear  functional  F,  namely, 

Rff,  fd  -r  V  IV  e-j<ra,t>  -  uK) 

B  (u>  v)  -  Jjd  2^KeK\2^(n,\)eBK  2  \E  ^(n.A)  _a;(K)| 

(v\e~i<n,t>  -t)*), 


|x(n,A)  —  x(K)\ 


■)  +  (C<“)fi„  v*)}dt 


(2.19) 


and 


F(v)  =  /  (f,  vK)dt. 

k£K  Ld 


Then  we  have  the  following  lemma. 


(2.20) 
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Lemma  2.6  Let  u,v  e  Hl{Q)  and  /  e  L2(Q),  and  let  u,v,  f  be  their  Fourier 
transform,  respectively.  Then  there  hold 

B{u,v)  =  (2TT)~dB{u,v)  (2.21a) 

and 

F(v)  =  (2ir)-dF{v).  (2.21  b) 

Proof.  For  v  =  Eme^  and  /  =  J2mezd  there  holds 

#(*)  =  E/  </..««)*  =  E  /  E 

kG/C  ld  kG/C 

—  (^0  ^  ^  ^  ^  ^  /m,K5  ^  • 

kG/C  rn£Zd 

This  leads  to  (2.21).  Similarly  we  have 

E  /  <c Wfl.,0.)*  =  (2»)dE  E  (Cw«m, «,*>■»,.}■  (2.22) 


It  is  easy  to  see  that 


f  W  E  I  /F(k,«,a)  (u\e~i<n,t>  -  uK)  (uAe  «».*>-(),) 

y/(1Z^  2^  2'  |a;(n,A)  _  ^(/c)  |  ’  |Jp(n,A)  _  x(*)|  ' 


zsG/C  (ti5A)6J5)c 


-  /  2^  2(E  2^ 

JI  K&K(n,X)eBK  m£Zd  1  1 

,kc _ A  7, 

Mn.A)  _  x(k)  I  ' 

m'eZd  1  1 

_  f  V  /p(K.n.A1  (um+n, A  ~  ^m.re)  (^m+rt,A  ~  vm,K,)~.  ,f 

12-^  2—J  2~i  2'  |x(n’A)  —  (c(K)  I  ’  |a;(n.A)  _ 

JI  KefC(n,\)€BKmeZd  1  11  1 

—  (9  V  V  V  ^  /pfc.n.A)  (Mm+n, A  ~  Miii,n)  (Vm+n, A  ~  um,/t)  \ 

—  v^71/  2  ->  2  j  2  ->  2'  |a.(m+n,A)  _  2;(m,/e)|  ’  |a.(m+n,A)  _  x(m,K)  I ' 


m£Zd  k&K  ( n,A)eB*  1  11  1 

which  together  with  (2.22)  yields  (2.21).  □ 

In  order  to  properly  define  a  variational  problem  over  Id  x  Kq, we  need  to 
introduce  new  function  spaces.  Let  L2(Kq)  and  Hl{, Kq )  be  the  spaces  of  grid 
functions  on  Kq  with  the  following  norms 

imIl2^)  =  ^2  i^i2 = ^2  ^2  i^i2 

«G/C  k£JC  KKs 


iHI HHKQ)  =  2^  2^  U(n,A)  _  g(«)|2  T-mL2(KQ). 

«GK(n,A)eBK  1  1 

L2  (ld,  H1(Kq))  and  L2  (/d,  L2(Kq))  are  spaces  furnished  with  the  norms  : 


-E  E 


wn,\  -  Wk 


+  im 


"L2(ld,W(KQ)) 


=  [  IN*)I 

Jld 


\2Hl(KQ)dt 
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and 

W^L^I^L^Kq))  =  fjd  IH*)II L*(KQ)dt- 
Lemma  2.7  If  /  G  Hl(Q),  then  f(t )  G  L 2  ( Id ,  H1(Kq )) ,  I  =  0, 1,  and 

II/WIIl*(«(M  =  (27r)<<ll/llflr'(e) 

Proof.  The  assertion  follows  easily  from  the  definition  of  the  spaces.  □ 

Remark  2.2  The  functions  u  =  u(t)  in  the  space  Hl(Id),  l  =  0, 1  are  vector  func¬ 
tions  with  sq  components,  and  the  functions  u  =  u(t,x M)  in  L2  (ld,  L2(Kq )) 
are  those  defined  on  a  semi-discrete  domain  Id  x  Kq  with  s  components.  Ob¬ 
viously,  the  space  L 2  (ld)  coincides  with  the  space  L2  ( Id,L2(KQ )),  and 

ll/llL2(/d,L2(ir<3))  =  ll/IU2(/d)- 

But  the  space  Hl{Id )  is  totally  different  from  the  space  L2  {ld,Hl(Ko)). 
The  latter  is  related  to  the  connectivity  BK,  and  the  former  is  not.  Further¬ 
more,  the  space  Hl(Id)  is  an  isomorphism  of  the  space  L\(Q),  and  the  space 
L2  (ld,  H1{Kq))  is  an  isomorphism  of  the  space  Hl{Q)  according  to  Lemma 
2.6. 

The  bilinear  form  B  in  (2.20)  and  linear  functional  F  in  (2.21)  are  defined 
on  L2  (ld,  H1(Kq))  x  L2  (ld,  H1(Kq))  and  L2  (ld,  H1{Kq)),  respectively.  The 
energy  space  E  =  E  (ld  x  Kq)  is  defined  as  one  equivalent  to  L2  (ld,  H1{Kq )) 
if  ^  0,  with  an  energy  norm 

\M\{I*XKQ)  =  B(w,wy'2. 

We  are  now  able  to  precisely  address  the  variational  problem  over  the  domain 
Id  X  Kq. 

Theorem  2.8  Let  B  and  F  be  the  bilinear  form  and  linear  functional  on 
L2  (ld,H\ Kq ))  x  L2  (ld,  Hl{KQ))  and  L2  (ld,H1(KQ)),  given  in  (2.19)  and 

(2.20),  respectively.  If  f  G  L2  (. Id,L2(KQ ))  and  C<*>  £  0,V«  G  K,  then  the 
variational  problem 

B  («,  v)  =  F(v),  V0  G  L2  (Id,  H\Kq ))  (2.24) 

has  a  unique  solution  0  G  L2  (/d,  H1(Kq)),  and 

ll“lli2(/d,ff1(^(?))  —  (L2(id,L2(KQ)))-  (2.25) 

Remark  2.3  Combining  Theorem  2.8  with  Theorem  2.3  and  Lemma  2.6,  we 
have  that  the  equivalence  between  the  problem  (2.10)  and  the  problem  (2.24), 
i.e.  the  equation  (2.24)  has  a  unique  solution  u  G  L2  (Jd,  ^(Kq))  and  the 
estimate  (2.25)  holds  for  /  G  (ld,  L2(Kq))  if  and  only  if  the  problem  (2.10)  has 
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a  unique  solution  u  G  Hl(Q )  and  the  estimate  (2.14)  holds  for  f  —  T  1(f)  G 
L2(g),  and  u  =  F-'iu). 

Applying  Fourier  transform  to  the  equilibrium  Equations  (2.6)  leads  to  an 
equilibrium  equations  over  Id  x  Kq 


_  ^  g(/c,n,A) 


(uAe  *<">*>  -  uK) 

[^(n.A)  _  3;(k)|2 


+  c  <*>««  =  fK,  Vk  g  /c. 


(2.26) 


Then  we  have  a  theorem  indicating  the  relation  between  the  solution  of  (2.24) 
and  the  solution  of  (2.26). 

Theorem  2.9  If  u  G  L2  (/d,  H1(Kq))  is  the  solution  of  the  variational  equa¬ 
tion  (2.24)  with  /  G  L2  (/d,  L2(Kq)),  then  it  satisfies  the  equilibrium  equa¬ 
tion  (2.26).  Vice  versa,  if  u  G  L2  (ld,  H1(Kq))  solves  the  equilibrium  equa¬ 
tion  (2.26)  with  /  G  L2  ( Id,L2(KQ )),  then  it  satisfies  the  variational  equation 
(2.24). 

Proof.  The  proof  is  analogous  to  that  for  Theorem  2.3.  □ 

2.4  Presentation  Formula  of  Solutions 
The  equation  (2.26)  gives  a  system  of  linear  equations 

<r{t)u(t)  =  f(t)  (2.27) 

where  u  =  (uf,  uj ...  Ug)T,  and  uK  =  (u\, . . .  usK)T ,  k  G  JC.  u  and  /  are  vectors 
with  sq  components,  and  <r  is  a  block  matrix 


a  ~  (<7’«»a)i</c,A<9  • 

Each  of  block  <rK\  is  a  s  x  s  matrix.  It  follows  from  (2.26)  that 

<7^=  (1  _e-i<n,t>)E(«^.«)  +  Cw,  (2.28o) 

ti£Bkk 


<tk,x  =  -  e-i<nd>E{K'n'x).  (2.28 b) 


Lemma  2.10  The  matrix  cr  has  following  properties  : 

(1)  <t  is  a  Hermitian  matrix  and  Hermitian  block  matrix; 

(2)  aK\(-t)  =  <rx,K(t)T  for  A,  /c  G  K,  and  <r(-t)  =  cr(t)T; 

(3)  a(t)  is  a  positive  definite  matrix  for  all  t  G  Id  If  ^  0  for  k  G  JC; 

.  (4)  If  =  0,  cr(t)  is  a  semi-positive  definite  matrix  for  all  t  G  Id  and  is 
a  positive  definite  matrix  for  all  t  G  Id  for  all  t  G  Id  except  t  =  0. 

Proof. 
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(l)It  follows  from  (2.28b)  and  (2.4)  -  (2.5)  that  for  k  ±  A 
(TK  x  =  —  e-i<n,t>g(/c,n,A)  _  _  E  e*<-n 

n€BK\  /  neB\iK 

—  _  ei<m,t>g(A,m,K)  _  — ^  ^ 

meBx,K  rneBx,K 


For  /c  =  A,  we  shall  write 
where 

B^K  ~  {P  =  (?li,  TI2  •  •  .  Tld)  G  Bkk  |  Tty  0, 1  ^  Z  ^  d} 

B KK  —  {n  =  (rii,  7i2  •  •  •  Tirf)  €  Bkk  |  Tty  ^  0, 1  5;  l  P  d} 

Due  to  the  definition  (2.3)  of  jBK)a,0  ^  BK>K,  which  implies  R+K  fl  R~>K  =j0.  By 
the  property  (2.4),  n  G  B^K  if  and  only  if  —n  G  B~K.  Therefore 

<rKfK  =  C(k)  +  4  ^  sin2  — E(K’n’K)  (2.29) 

orKiK  is  a  real  matrix.  Thus  we  have  shown  that  a  is  Hermitian  block  matrix. 
Note  that  C M  and  E^K,n,A^  are  symmetric  matrices,  which  implies  the  cr  is  a 
Hermitian  matrix  as  well. 


(2)  For  A  7^  k,  by  the  properties  (2.4)  and 


o’k.aH) 


-  5Zn6BK,  a 

n£B\iK 


e~i<m,t>  g  A) 


J2-n€Bx,K 

~  ^ meBx,K 


e-i<-n,*>  jj(/c,n,  A) 


It  is  trivial  by  (2.29)  that  erK>K(t)  =  <rK)/c(— t),  Vk  G  /C.  Since  each  block  <r\,K{t) 
is  symmetric,  we  have  er(— t)  =  <r(t)T. 

(3)  Let  6(t)  =  (fc(t))  G  L2(Id,  Hx{Kq),  and  let  b  G  P~l  (&(t)).  Then,  by 
Lemma  2.6,  b  G  ff 1(Q),  and 


{a{t)b(t),b(t))dt  =  B(b,b) 


=  E  E  E  <E<“,”'A)  (*»+%»  -  W)  .  (f'm+n.A  -  W)> 

meZd  k£Kq  (n,A)€B« 

+  y^(C^6m+w, bm+K)  >  di\b\2Hi^  +  d2||5||i2(g) 
where  d\  and  d2  are  given  in  (2.13).  Note  that 

=  [  \Ht)\w(KQ)dt>  IHIi,2(a)  =  f  \\Ht)\\L2(KQ)dt’ 

J id  J  ld 

which  implies  that 

J^(<r(t)b(t),b(t))dt  >  d\  \kt)\m(KQ)dt  +  d2  WKtfWh^dt- 
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and  that  for  almost  every  t  G  Id 

{a(t)b(t),b(t))  >  +  d2\\b{t)\\2L^Kq)  >  d2\\b(t)\\h(KQ)- 

with  c?2  >  0.  Note  that  cr(t)  is  a  Hermitian  matrix  and  analytic  in  t.  Hence 
<j{t)  is  positive  definite  for  all  t  G  Id. 

If  =  0  for  all  k  €  1C,  there  holds 

(<r{t)b{t),b(t))  >  di\b(t)\2HHKQ)  >  0. 

cr(t)  is  semi-positive  definite  for  all  t  G  Id.  By  [?],  det(<r(t))  =  0  if  and  only  if 
t  =  0,  which  implies  that  a(t)  is  positive  definite  if  t  ±  0.  □ 

In  terms  of  <r_1(t),we  now  can  derive  a  representation  formula  of  the  solution 
the  lattice  problem. 


Theorem  2.11  If  ^  0  for  k  €  K,  and  /  €  L2(Q),  the  solution  of  the 
lattice  problem  (2.6)  can  be  represented  by 


u  =  T  1  (<7  x(£)/(£)) 

(2.30a) 

where  f  =  F(f)  E  L2(Id),  with 

um  =  (27r)_d  f  (T~l{t)f(t)e~l<m,r>dt 

Jid 

(2.306) 

and 

ll^-1  (^W-1/^))  II H'{Q)  <  C\\f\\L2^d>L2^KQ)y 

Proof.  Since  cr(t)  is  positive  definite  for  all  t  G  Id  if  ^  0, 

u  = 

is  the  solution  of  the  equation  (2.26).  By  Theorem  2.8,  it  solves  the  variational 
equation  (2.24),  and 

INI  L2(ld,m{KQ))  —  (^II/IIl2(/<<,l2(Kq))- 
Let  u  =  T~x  (• u(t )).  Then  u  G  Hl(Q)  and  solves  the  equation  (2.6),  and 


\\T  1  (<r(i)  1f{t)'j  ||hi(c?)  —  IMIi/qe)  ^  ^INI L2{id,Hi(KQ)) 

-  c\\f\\L2(idtL2(KQ)y 


□ 


If  =  0  for  all  k  G  1C,  then  for  any  b  G  E(Kq )  Since  <r(t)  is  a  positive 
definite  matrix  for  t  G  Id  and  £  /  0,  <r_1(t)  exists  for  all  £  /  0  and 

u(t)  =  <r_1  (£)/(£). 

i T~l{t )  is  singular  at  the  origin,  and  u(t)  has  pole  at  t  =  0  in  general  if  f(t)  G 
L2(Id)( i.e.  /  G  L2(Q).  Therefore,  the  integral  in  (2.30)  may  diverge  due  to  the 
singularity  at  the  origin.  To  make  the  integral  converge,  we  have  to  impose 
additional  conditions  on  /(£ )  which  are  able  to  absorb  the  singualrity  of  a  at 
the  origin.  For  instnace,  /  G  Hper  (ld)  with  u  >  1  and  vanishing  at  t  =  0, 
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which  is  equivalent  to  /  €  L2(Q)  with  v  >  1  and  Yhmqzd  ^2k<ek  /»,«  =  0-  Based 
on  the  analysis  on  the  lattice  probloems  without  absolute  terms  in  next  section, 
the  following  theorem  on  the  existence  and  uniquness  solution  and  validation 
of  the  representation  formula  can  be  proved. 

Theorem  2.12  If  f  e  L2{Q)  with  u  >  1  for  d  —  2  and  v  —  1  for  d  =  1, 3, 
and  Y!m<zZd  S/teJC  /n,K  =  0,  then  the  lattice  problem  (2.6)  with  =  0  has 
solution,  which  can  be  representd  by  the  formula  (2.30)  and 

IMUms)  =  \\f~1  (^W-1/^))  I \hhg)  <  C'll/ll^er(/d)  —  WfWmg)- 


2.5  Two  Examples  Of  Lattice  Problems 

We  will  analyze  two  concrete  lattice  problems.  One  is  one-dimensional,  and 
another  is  two-dimensional.  Although  the  structures  of  these  two  problems 
are  simple,  the  analysis  we  carry  out  here  can  be  generalized  to  other  lattice 
problems. 

A  lattice  problem  in  one  dimension 

Suppose  elastic  rode  of  two  different  materials  with  half-unit  length  and  in¬ 
tersection  area  A  are  connected  by  hinges  at  nodes,  shown  in  Fig.  2.1.  The 
master  cell  Q  =  [0,1),  containing  two  nodes  x^k\k  €  /C  =  {1,2}.  The  nodes 

in  cells  Qm  for  m  €  Z  are  denoted  by  x^m,K\  K  =1,2.  Let  M  —  {x^m,K\m  € 

Z,k  =  1, 2}  denote  the  global  mesh  containing  all  nodes.  Suppose  that  the 
rods  are  furnished  with  springs  at  each  node.  By  E\  and  Ei  we  denote  the 
Young’s  modulus  of  the  rode,  and  by  A.C\  and  4C2,  the  Hook’s  constant  of  the 
springs,  respectively.  A  lattice  Q  denotes  such  a  structure,  connectivity  and 
periodic  translation. 

Equilibrium  Equation 

Let  Uj  and  4  fj  denotes  the  displacement  of  the  rods  and  external  force  at  the 
nodes  x we  have  the  following  equilibrium  equations 

A{E\(um> 2  Um,i)  +  #2 Km  1,2)}  T  =  fm,  1  gjN 

— A{E2(um+iti  —  Umfi)  +  Ei(um>2  —  Mm,i)}  +  C2«m,2  =  fm,  2 


Variational  Equation 

The  corresponding  variational  equation 

B(u,  v )  =  F(v). 

where  u,  v  and  /  are  functions  defined  on  Q,  with  the  bilinear  form 

B(u ,  u)  =  y^^g?(Um.2  Wm,x)  AE\  (um,2  T 

T  Umfi)  AE2{vm+iti  Vm,2)  +  £2^777,2^771,2 

and  the  linear  functional 

F{v)  =  fm,jvm,j 

meZ  j=l,2 
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The  energy  space  E(Q)  contains  functions  on  Q  with  finite  energy  E(u), 

E( u )  =  u)  =  |  YlmeZ  AEi(um,2  —  wm,l)2  +  CWm.l 

“I”  2  5- im€Z  AE2{um+l,\  ^m,2)  ”h  C2Um2- 

The  spaces  Hl(Q)  and  L2(Q)  are  furnished  with  the  norms 

IMIffi(£)  =  Mlf1(<?)-HMIi2(0)  =  y~l  (um, 2  ~  Um,l)2  +  (Urn+1,1  ~  ttm,2)2  +  |Mll2(g) 

m£Z 

and 

iwiiw  =  E  + k*i2- 

TYlCZ 

Hence,  the  energy  norm  ||u||e(0)  =  E{u)1!2  is  equivalent  to  the  norm  of  the 
space  Hl(Q)  if  C\  +  C2  /  0,  and  equivalent  to  the  semi-norm  of  the  space 

H\G)  ifC!  +  C2  =  0. 

Fourier  transform 

For  /  =  f(m),m  e  2,  we  introduce  the  Fourier  transform 

E(f)  =  f(t)  =  ^2f(m)eimt,te  I  =  (-7T,  tt) 
n 

which  realizes  an  isomorphism  between  L2{Q)  and  L2(7),  and  between  L2{Q) 
and  H"I),  where  the  space  L2(Q)  is  defined  as  a  weighted  space  with  a  weighted 
L2-norm 

IMIl2(G)  =:  ^  +  m2)U(\Um,l\2  +  lMm,2p)- 

m€Z 

The  inverse  Fourier  transform  gives  /  =  E~x(f)  with 

Applying  the  Fourier  transform  to  the  equations  (2.31),  we  obtain 


A(Ei  +  E2)ui  —  A(E\  +  E2eu)u2  +  C’i'Ui 
— A{E\  +  Eke-*)*!  +  A{Ei  E2)u2  +  C2u2 
The  corresponding  matrix 


<T 


En 

—E2i 


—E12  \ 
E22  ) 


(2.32) 


with  E\\  —  A{E\  +  E2)  +  Ci,  E22  —  A(Ei  +  E2)  +  C2,  E\2  —  A{E\  +  £l2e*4)  and 
E2i  —  A(Ei  +  E2e~tt).  <r(t)  is  a  Hermit  matrix,  and 

det(a )  =  £A? EiE2sin2t/2  +  A{Ci  4-  C72)(i?i  -1-  E2)  +  CiC2. 

Obviously,  if  C\  +  C2  >  0,  a(t)  is  positive  definite  for  tel,  cr~l{t)  exists, 


<T  x(  t ) 


1  /  E02 

det((T )  \  E12 


E2 1  \ 
E01  ) 


Representation  Formula 
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If  Ci  +  C2  >  0,  the  solution  of  equation  (2.32)  can  be  given  by 

( ^ )  =  (  a  )  • 

Therefore,  the  solution  of  the  problem  (2.31)  can  be  represented  as 

«.  =  *  =  1.2 

1<^<2 

with  x 

^1,1  (i)  =  detect)))-®102 ’  02,2  fa)  =  rfet^r(t)))'®'01’ 

01>2fa)  =  det(<T(t))^21>  02, lfa)  =  dei(<r(i))-®'11- 

If  Ci  =  C2  =  0,  <rfa)  is  positive  definite  for  t  €  /  except  t  =  0,  and  has  a 
pole  of  order  2  at  t  =  0.  If  |/(t)|  =  0(|t|a)  with  a:  >  1/2  near  the  origin,  then 
u  has  a  pole  of  order  1  at  t  =  0.  Hence  u  £  L2(I ,  L2(Kq)),  but 

(<r(t)u,u)  =  (<r_1fa)/,/)  <  00 

which  implies  that  |M|e(£)  <  00  if  /  €  L\(Q)  and  Ymez  XLg/c  /m,«  =  0. 
Hence  the  problem  (2.31)  with  Ci  =  C2  =  O'  has  a  unique  solution  in  E(Q)  up 
to  a  constant.  The  argument  above  based  on  Fourier  transform  can  be  utilized 
for  periodic  lattice  problems  in  entire  spaces  without  absolute  terms  [?,  15], 
and  leads  to  a  result  coinciding  with  Theorem  2.16. 


A  lattice  problem  in  two  dimensions 

Suppose  elastic  rods  with  intersection  area  A  and  unit  length  are  connected 
by  hinges  at  nodes  xkj  —  ( k,j )  €  Z2,  periodically,  shown  in  Fig.  2.2.  The 
master  cell  Q  =  [0,  l)2,  in  which  there  is  only  one  node  x^\  and  the  index  set 
K  =  {1}.  The  mash  M  =  U W=(*j)esa*(m,1)  =  U(kj)ez*Xkj  =  U (k,j)ez*(k,j). 
We  will  use  Xk,j  to  denote  the  nodes  in  stead  of  x^m,1\  Suppose  that  the 
rods  are  furnished  with  springs  at  each  node.  By  E  and  C  =  diag(C,C), 
we  denote  the  Young’s  modulus  of  the  rods,  and  the  Hook’s  constant  of  the 
springs,  respectively.  A  lattice  Q  in  two  dimensions  denotes  such  a  structure, 
connectivity  and  periodicity. 


Equilibrium  Equation 

Let  ukj  =  fa4j,4j)  and  fk,j  =  ( )  be  the  displacement  vector  and  ex¬ 
ternal  force  vector  at  the  node  xkj  =  (k,j).  We  have  the  following  Equilibrium 
equation 


(  +  ^12Uklij_i) 

(-A12  ujj  +  AtttijJjj.j) 


(2.33) 


where 


A1  =  4+1, 3  -  uk!p  Aatig  =  1  -  ui%  A12 =  (4+1J+1  -  n^/V2. 


Variational  Equation 
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The  corresponding  variational  equation  is 

B(u,v)  =  F{v ) 

where  the  linear  functional  F  and  the  bilinear  form  B  are  defined  as  ' 

F{v)=  fliVkj 
(k,j)ez2 

and 

B(u, v)  —  +  (&2uk,j) 

+|(A12Ufcj)TE*A12Mfcj  +  Cv%jVkj 

where  E  =  AEX  and  E*  =  AEbbT,  I  is  an  identity  2x2  matrix,  and  b  is  a 
vector  =  (cos|,sm|)T. 

Fourier  Transform 

We  introduce  the  Fourier  transform  for  functions  /  —  {/*j,  (k,j)  e  Z2} 

m = m  =  E  a ,jei(“,+>fe),  t = (tu  *2)  e  (-Jr, 

(k,j)€Z2 

which  leads  to  an  equation  in  matrix  form 

er(t)u{t)  =  f(t )  (2.34) 

where 

<r(t)=(aa  Cl2) 

\  U21  <722  ) 

with 

(7n  =  AE(4sin2^  +  y/^sin2^^-)  +  C 
a 12  =  <721  =  — \/2  AEsin2^}^ 

a22  =  AF(4sm2f  +  \/2sm2^2)  +  C 

<r(t)  is  a  real  and  symmetric  matrix,  and 

det(<r )  =  T2i?2(16sm2^-sm2^  +  4\/2  sin2 — ^  (sin2 1-  +  sin 2^)) 
+2AFC(2sm2|  +  2sm2|  +  v^sm2^)  +  C2. 

If  C  >  0,  the  matrix  cr-1^)  is  analytic  in  a  strip  E«5  =  {2  :  |/m2|  <  5}  with 
<5  >  0,  and  <r(t)  is  positive  definite  for  t  €  (— 7T,  7 r)2,  and 

^  _<T21  )  • 
v  '  det{<r)  \  —<7\2  0u  J 

Representation  Formula 

If  C  >  0,  then  cr(t)  is  positive  definite,  then  we  a  solution  of  the  equation 

(2-34) 

(ii(1)((),'S<2,(())r  =  o’  1(*)(/<1>(i)./(2>(*))T- 

Then  solution  of  the  equation  (2.33)  can  be  represented  as 
(u<i),u(2))t  =  ?-\<t 

=  ^-l(E.<2<2  (*).  E,<2<2  WF- 
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with 

*,.(<)  =s^-)(AE(i3iv?ti  +  V23mi^)+C), 

M*)  =  +2,1  (*)  =  dfrT/UEmttopL, 

■M*)  =  +  \/2Sm2^)  +  C). 

If  C  =  0,  cr(t)  is  positive  definite  for  t  €  I2  except  t  =  0  and  has  a  pole 
of  order  2  at  t  =  0.  If  |/(£)|  =  0(\t\a)  with  a  >  0  near  the  origin,  then  the 
equation  (2.34)  has  a  solution  u(t)  £  L2(I,L2(Kq)),  and  has  a  pole  of  order  1 
at  t  =  0.  It  can  be  shown  that 

(cr(t)u(t),u(t))  =  (a--1  (t)f(t),f(t))  <  oo 

which  implies  that  IMIe(0)  <  oo  if  /  €  L\{Q)  and  J2meZ2  jc  /»»>«  =  0  and 
that  the  problem  (2.33)  with  (7  =  0  has  a  unique  solution  in  E{Q)  up  to  a 
constant.  This  coincides  with  Theorem  2.16 

3.  UNSTRUCTURED  LATTICES  IN  ENTIRE  SPACES 


3.1  Non-periodic  lattices  in  entire  spaces 

A  lattices  Q  in  R d,d  =  1, 2, 3  is  called  unstructured  if  it  is  not  periodic  or 
quasi-periodic.  It  is  comprised  of  countable  number  of  nodes  distributed  in 
entire  spaces,  which  are  connected  by  elastic  rods.  Let  K  =  {xk,  k  G  Af  = 
{1,2.  . . . }}  denote  the  set  of  nodes,  and  let  be  elastic  rod  connecting  the 
nodes  Xk  and  xi  with  the  Young’s  modules  and  intersect  area  A.  The 
distribution  and  connection  of  the  nodes  can  be  periodic  or  non-periodic.  The 
assumptions  on  the  connectivity  and  rigidity  for  periodic  lattices  in  previous 
section  should  remain  the  same,  but  given  in  the  non-periodic  setting.  The 
assumption  on  connectivity  is  now  decribed  as  follows: 

(CM)  Each  node  is  connected  to  others  by  the  rods, at  least  one  node  and  at 
most  M  nodes. 

(C".2)  The  length  of  rods  are  uniformly  bounded  from  below  by  b\  and  from 
above  by  62,  i-e.  for  any  Xk  and  Xe  which  are  connected,  there  holds 

h  <  \xk  -xt\  <  b2. 

(C1. 3)  Any  two  nodes  Xk  and  xg,  k,£  €  Af  are  linked  by  a  chain  Lk/.  x[k'^  = 
Xk  — >  x^'^  — y  xf’^  —>■•••—>•  x^’^  =  xg  such  that  x[k’^  is  connected  to 
x[+f  \  1  <  t  <  s  —  1,  and 

I xk  -Xi\<  ^2  \xt+i  ~  xtk'e) I  <  V\xk  ~  xt\. 

KKs-l 

There  may  be  several  chains  connecting  the  nodes  Xk  and  xg,  Lk,g  denotes 
always  the  shortest  one.  In  particular,  s  =  1  if  Xk  and  xg  are  connected. 


To  describe  effectively  the  connection  of  the  nodes,  we  introduce  a  set  of 
indices 


Kk  =  {£  £  M\ xg  is  connected  to  Xk} 
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which  together  with  the  assumption  (C.1)-(C.3)  gives  the  connectivity  of  a 
lattice.  The  nodes  {a;*,  k  €  Af}  and  the  connectivity  of  the  nodes  uniquely 
characterize  the  structure  of  a  lattice  denoted  by  Q  =  G{K,  K.k). 

For  sake  of  convenience,  we  partition  the  entire  space  Rd  into  cells  Qm,  m  = 
G  Zd  such  that  they  are  disjointed  and  U m^zdQm  =  Rd. 
Here  Z  —  {0,  ±1,  ±2 . . . }  as  in  Section  2.  Qm  is  an  interval  in  R1,  a  polygon 
in  R2,  and  a  polyhedron  in  R3.  Qm  D  Qn  is  avertex,  or  an  edge,  or  a  face, 
or  empty  set  for  any  n,m  G  Af,n  ±  m,  where  Qm  denotes  the  closure  of  Qm. 
Furthermore,  a  partition  of  Rd  is  assumed  to  satisfy 

(Q.l)  Vertices  of  Qm  are  nodes  of  Q,  and  edges  of  Qm  are  elastic  rods; 

(Q.2)  Every  node  must  be  a  vertex  of  a  cell,  or  is  located  on  an  edge  of  a  cell, 
or  in  interior  of  a  cell,  and  each  node  is  shared  by  at  most  Uq  cells; 

(Q.3)  The  diameter  of  cells  is  uniformly  bounded, 

d q  ^  diag(Qm)  <  <1q] 

(Q.4)  The  number  qm  of  nodes  in  each  Qm  is  uniformly  bounded  by  q ; 

(Q.5)  Each  node  in  a  cell  is  connected  to  at  least  one  node  in  the  same  cell; 

(Q.6)  Each  cell  is  rigid.  According  to  [1],  there  is  a  sufficient  condition  for  the 
rigidity  nof  a  cell, 

Cfn  dqm  d(jl  “I- 1  )/2 

where  em  be  the  number  of  rods  in  the  cell  Qm,  which  leads  to  rigidity  of 
lattices. 

We  may  introduce  a  local  number  (m,  k),  k  €  Afm  =  {1)2,...,  qm}  for  the 
node  Xk  in  the  cell  Qm  such  that  Xk  =  x^m,K^  and  \xk\2  =  \m\ 2  =  Yli<i<d  \mi\2- 
The  connection  of  nodes  can  described  locally  described  by 

B™  =  {(n,  A)  |  x(n+m>A)  is  connected  to 


Proposition  3.1.  If  Xk  —  x^m,K^  and  xi  =  x('m,x^  are  in  the  same  cell  Qm,  the 
number  of  nodes  on  the  chain  Lki  is  not  more  than 

0\ 


Proof.  Due  to  (C'.3),  there  is  a  chain  Lk,i  linking  the  nodes  Xk  and  Xg  :  x\  g  = 
Xk  xle  — >  x\g  xsk>g  =  xg.  By  (C.2),(C.3),and  (Q.3),  there  holds 


sh  <  ^2  1 4+?  ~  xtk’e) I  <  V\xk  -xi\<  vdQ, 

l<t<s-l 


which  leads  to  the  assertion  of  the  proposition. 


□ 


Let  Uk  =  (itjfe,  . . . ,  usk)r  be  a  vector  function  on  the  nodes  {xk,  k  G  A f}.  As 
in  previous  section,  we  shall  deal  with  the  problems  with  s  =  lforl<d<3 
and  the  truss  problems  with  non-rigid  joint  in  one,  two  and  three  dimensions 
for  which  s  =  d  =  1, 2, 3.  We  furnish  the  rods  with  springs  in  the  axis  directions 
at  each  node  with  Hook’s  coefficients  denoted  by  diagonal  matrices  Ck,k  e  Af. 
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We  assume  that  the  ratio  of  the  length  of  the  rods  and  the  intersect  area  A  of 
rods  >>  1. 

If  external  forces  exert  on  the  rods  at  the  nodes  xk,  denoted  by  /  =  {/*,,  k  e 
A/"},  we  have  an  equilibrium  equation  for  the  elastic  rod  problem 


_  £  E(*J)  (U‘  u*l  +  ckuk  =  /fc,VAt  G  M 


t€Kk 


(3.1) 


where 


and 


e(M)  _  AE(k,t)  £*)  {Ei 

l^-Xfcl2  |rPf-a:fc|2 


Cfc  =  dia<?(Cf,  C2fc, . . . ,  Ckd),  Cf  >  0, 1  <  l  <  d. 

Let  Hl(Q)  and  L2{Q)  be  the  discrete  Sobolev  spaces  with  the  norms 


and 


imiW)  =  E  w2 

Ml/rqc)  ==  M  h1(6)  +  IMIi2(0) 


where  \u\Hi(g)  is  the  semi-norm, 


\u\hHG)  ~  XI 


\ui  ~^fc|2 
\xe-xk\2' 


(3.2) 

(3.3a) 

(3.36) 

(3.4) 

(3.5a) 


k&MieKk 

Similarly,  the  corresponding  variational  equation  can  be  derived 

B(u,v)  =  F(v)Vv  e  H\g) 

where 

is  a  bilinear  form  on  Hl{Q)  x  Hl(Q),  and 

fee  Af 

is  a  linear  functional  on  Hl(Q). 

The  energy  of  the  lattice  Q  is  defined  as 

G(„)  =  B(u,  u)  =  E(  E  +  <C V  *»  (3.6) 

The  energy  space  denoted  by  E(Q)  is  the  family  of  all  grid  functions  u  on  Q 
with  finite  energy  G(u)  defined  in  (3.6),  and  G(u)1!2  is  referred  as  the  energy 
norm  ||w||b(^).  The  assumption  (Q.6)  leads  to  the  rigidity  of  lattices,  and 

EJ>«^S>=o- 

if  and  only  if  u  is  a  rigid  body  motion. 
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At  the  node  xk  =  x^m,K\  we  may  write  uk  =  umjK,Ck  =  C (m,K\fk  =  fm 
etc.  Due  to  the  disjointedness  of  cells,  we  have  the  following  proposition  which 
allows  us  to  write  the  norms  either  in  global  numbering  or  in  local  numbering. 


Proposition  3.2  For  u  E  L?(G),  v  E  Hl(Q )  and  w  E  E(@),  there  hold, 
respectively, 

IMIl2(£)  =  \Um,K^  1 

meM'  K€A/*m 

iinii2  =  v  I^.a  z  ^l2 

"  ^H1(Q)  /  j  ^  /  j  |^(n,A)  _  /p(m,/c)  |2  ’ 

meN  KZNm  {n,\)eK™  '  1 

and 


_  1  /T ^ n ^  ^n}A  v 

2  2-^meAf  2-«tK£Afm  ^-'(n,A)e/C™  \  |^(n,A)  _  /p(m,/c)| 5  |^(n,A)  _  %(m,K)  | ' 


+  Erne*  £k6M»  (CMV,  Wm,*>  • 

As  analogue  of  the  Theorem  2.3,  we  have  the  existenc  and  uniqueness  for 
the  non-structured  lattice  problems  (3.1). 

Theorem  3.3  The  variational  form  B  on  Hl{Q)  x  Hl(Q)  given  in  (3.6)  is 
continuous  and  coercive  if  Ck  ^  0  in  the  sense  that  there  exists  a  node  xkm  = 
x(m,Km)  jn  eacjj  cen  such  that  Ckm  ^  0,  and  the  problem  (3.1)  has  a  unique 
solution  u  €  Hl{Q). 

Proof  The  proof  is  analogous  to  that  for  periodic  lattice  problem  (2.6).  For 
the  detail  in  the  non-periodic  setting,  we  refer  to  [?]. 

3.3  Problems  Without  Absolute  Terms 

The  proof  for  the  existence  and  uniqueness  of  solutions  for  non-periodic 
lattices  with  absolute  terms  is  similar  to  those  for  periodic  lattices  which  we 
have  invetigated  intensively.  We  should  concentrate  on  the  problem  of  non¬ 
periodic  lattices  without  absolute  terms  in  this  section  for  which  the  Fourier 
transform  can  not  be  used.  Fourier  transform  is  an  extremely  powerful  tool  for 
periodic  lattices  with  and  without  absolute  terms.  It  is  worth  to  pointing  out 
that  in  practical  applications,  problems  of  structured  or  unstructured  lattices 
are  associated  with  no  absolute  term,  i.e.  =  0  for  all  k  €  M.  We  need  to 
develop  a  new  approach  for  analysis  of  non-periodic  lattices  without  absolute 
terms, which  is  extremely  important  practically  and  theoretically 
We  seek  u  E  E{Q)  such  that 

B(u,v)  =  F(v)\/vEE(Q)  (3.8) 

where 


E(v)  =  (f,v)g  =  ^(/fc,Ufe). 

keN 


and 


(3.96) 
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For  the  problem  (3.8),  the  energy  space  E(Q)  does  not  coincide  the  Hl{Q), 
and  E{Q)  is  not  embedded  in  L2{Q).  Actually,  the  “energy  norm”  ||u||£(0)  is 
equivalent  to  the  semi-norm  \u\h^(q)-  Consequently,  the  solution  of  the  lattice 
problem  is  not  unique  and  may  not  exist  in  Hl(Q).  Hence  we  have  to  modify 
the  space  Hl{Q ),  such  that  the  solution  exists  uniquely  in  modified  space 
Hl(G)  with  the  norm  equivalent  to  semi-norm  of  Hl(Q). 

The  modification  of  the  spaces  and  the  equivalence  of  the  norms  for  one 
dimension  is  different  from  those  for  two  and  three  dimensions.  We  shall 
address  them  separately.  To  this  end,  we  define  a  weighted  space  Lla(G)  for 
all  dimensions  with  the  norm 

INI ha{g)  =EfeeA7(1  +  l^lTlog2ff(l  +  M  kl2  . 

=  Em€sdEK€A7ji  +  HTiog2CT(i  +  H)K«|2  ^  ; 

where  v  and  a  are  real  numbers.  We  shall  write  Ll0(Q)  =  L2(Q),  L2v0{G)  = 
Ll(G).  Obviously,  Ll^Q)  D  L\Q)  iiv,a  <  0,  and  C  L\Q)  if  i/>  >  0. 

3.3.1  Problems  without  absolute  terms  in  one  dimension 

Lemma  3.4  If  v  G  E(G),  and  v0  =  0,  then 

IMiLijta)  <  C\v\Hi(g).  (3.11) 


Proof.  Due  to  the  proposition  3.2,  we  shall  use  local  numbering  for  the  nodes. 
First,  suppose  that  and  in  the  cell  Qm  are  always  connected  for 

any  k  G  Mm  with  k  ^  1.  Then 


Emez  E 


umiK\ 


TtlE.2  ^  ^2 


—  ^  [  ^rrneZ  ^m,l  )2  +  E 


l«wi 


meZ  1  +  m2 


(3.12) 


If  there  are  some  k  G  A fm  such  that  x^m,K^  and  x^m'1^  are  not  connected, 
there  always  exists  by  (C'.3)  a  chain  :  x^m,r>  —  rr(ni,Al)  -4  x^712^21  •••—»• 

x(ns,\s)  —  and  due  to  Proposition  3.1 


£  K«  -  <wl2  <  c  £  £  £  <  ciHiptc- 

rnez  m.ez  KEArm  (n,A)eBJ? 

(3.13) 

Hence,  (3.12)  holds  for  the  cases  that  x^m,K^  and  are  connected  or  not 
connected. 

Let  wm  =  vm^.  It  is  sufficient  to  show  that 


E 


\Wm£ 
1  +  m2 


00 

^  C  ^  \wm  u>m-i| 

m= 1 


2 


(3.14) 


We  assume  that  v0  =  wo.i  =  w0  =  0.  We  have  by  Cauchy  inequality 

oo  |  i2  00  1  /  m  \  2 

ElT^  S  Einb  EH-h) 

m=  1  m=l  \  j=l  / 

°°  i  /  m  \  /  m  \ 

£  ErTrffE^-^nfEr'J 

m=l  \J=1  /  \J-1  / 


where  £  G  (0, 1),  arbitrary.  Note  that 


m  pm 

]Tre<c'/  cedt;<cm l- 

.7=1 


which  implies  that 


oo  i  m 

nr  e 


<  c'^K  -^-i)2ieziy 

j=l  m=j 


For  j  <  1,  there  holds 


OO  _£  OO  /»oO 

Ei r~«< or 


which  leads  to  (3.14).  Therefore,  there  holds 


oo  2 

V^  Vmi 


Similarly,  it  is  true  for  m  <  0.  Therefore  we  have 


y 

1 1  +  m2  1 

m(zZ  m€Z 


(3.15) 


Note  that  um)i  and  um-i,i  may  be  connected  or  not  connected,  the  argument 
for  (3.13)  can  be  carried  once  more  here.  Hence,  we  have 

~y  ]  |^m,l  ~  vm- 1,1 1  ^  ^  y  }  ]  ^  1  |^m+n,A  ^m,re|  <  @  » 

me.Z  mez  KeMm  (n,\)€B™ 

which  together  with  (3.12)  and  (3.15)  leads  to  (3.11).  □ 

Theorem  3.5  If  /  €  L\{Q)  and  //ft  =  0>  then  for  any  v  €  E(Q)  it  holds 
that 
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Proof.  Note  that  Ll(Q)  c  L\(Q)  and  the  sum  YlkeM  fk  exists.  By  Lemma 
3.4  we  have 


< 

< 


\(Mo\  = 


£  fkVk 


keAT 


£  /*  0 V k  -  Vo) 


I  feejV 
\  1/2 


£(i  +  W2)  IM2  £  T 

ueV  /  ' 

c'll/llLf(e)lwlH1(e) 


(v*  -  «b) 


2  \  1/2 


^keAf 


+  M 


□ 


We  introduce  the  space  Hl{Q )  with  the  norm 

IMIhhs)  =  {iMIli^o +  > 

and  a  quotient  space 

H\Q)  =  H\G)/V r 

with  norm 

hWano)  =  ™ij\u  -  a\\ff1(6) 

where  Vr  is  the  set  of  all  constant  functions  on  Q,  which  is  a  subspace  of  Hl(G). 
Then,  by  Lemma  3.4  , 

IMIifi(a)  —  \u\hhg)  -  IMU(s)- 

In  the  framework  of  the  space  Hl{G)  and  Hl{G)  we  are  addressing  the 
existence  and  uniqueness  of  the  solution  of  the  problem  (3.8)  with  d  —  1. 

Theorem  3.6  If  /  €  L\{G)  and  Y^kzN  fk  =  0>  then  the  problem  (3.8)  with 
d  =  1  has  a  solution  u  €  E(Q),  and 

M HHQ)  <  C\\f\\L{(g)-  (3-17) 

The  solution  is  unique  up  to  a  constant. 

Proof.  Due  to  the  equivalence  between  |w|#i(£)  and  IMIjyi(g) 

\B(u,v)\  <  C\u\Hi{g)\v\HHg) 

—  C|Mljfi(e)IMIff(e) 

and 

\B(u,u)\  >  D\u\2Hi^  >  D\\u\\fyi(gy 
By  Theorem  3.5,  it  holds  that 

By  Lax-Milgram  Theorem,  the  variational  problem  has  a  unique  solution  u  € 
H\g),  and 

which  implies  (3.17)  and  the  uniqueness  of  the  solution  in  E(Q)  up  to  a  con¬ 
stant. 
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3.3.2  Problems  without  absolute  terms  in  two  dimensions 

With  suitable  spaces  L2pa(Q)  and  Hl{Q)  we  are  able  to  address  properly  the 
problem  (3.8)  without  absolute  term  in  two  dimensions.  The  next  theorem  is 
essential  to  the  existence  and  uniqueness  of  solutions. 

Theorem  3.7  If  /  <E  L\i{Q),  and  Y^ketf  =  tlien  for  any  v  e 

l(/»c|  <  C\\f\\L>js)\v\HHg)  (3.18) 

with  constant  C  independent  of  /  and  v. 

The  theorem  is  parallel  to  Theorem  3.6  for  one  dimension,  but  the  proof  for 
two  dimensions  needs  some  embedding  results  for  functions  in  R2,  which  is 
contained  in  Appendix  A.  In  order  to  use  these  results,  we  have  to  extend  by 
linear  interpolation  grid  functions  on  Q  to  whole  space  R2. 

Let  Km  be  a  set  of  nodes  which  are  located  in  Qm  where  Qm  is  the  closure 
of  the  cell  Qm.  Obviously,  x^m,K^  E  Km  for  all  k  E  1C,  and  some  nodes  a^^  in 
neighboring  cells  are  included  as  well.  Let  K^,  and  be  subsets  of  Km 
for  nodes  at  vertices,  on  edges  (not  including  vertices)  and  in  the  interior  of 
Qm,  respectively.  Then  Km  =  U  U  K!m. 

By  Tm  =  {ti,  i  =  1, 2, . . . ,  T}  ,  we  denote  a  triangular  partition  of  Qm 
satisfying  the  following  conditions: 

(T.I)Lt^  =  Km,  where  Vjm  denote  a  set  of  all  vertices  of  the  partition  Tm; 

(T.2)  The  partition  is  regular,  i.e.  U  n  tj  for  j  is  a  vertex,  or  a  whole  edge, 
or  empty. 


Fig.  3.1  Triangular  partition  of  cell  Qm 


The  construction  of  such  a  partition  can  be  started  with  an  initial  partition 
7^  of  Qm  for  which  (T.2)  holds  and  Vjv  =  K^n,  shown  in  Fig.  3.1(a).  For  the 
partition  7^  there  may  be  some  nodes  in  the  interior  of  triangles  tjS.  If  a  node 
x(n,K)  g  js  jn  the  interior  of  ti,  we  divide  ti  into  three  smaller  triangles  by 
connecting  £(”>*)  to  four  vertices  of  ti,  shown  in  Fig.  3.1(b).  Repeating  the 
process  for  each  node  in  the  interiors  of  all  simplices,  we  have  a  partition  7^ 
of  Qm  for  which  (T.2)  holds  and  no  node  G  Km  is  located  in  the  interiors 
of  all  triangles.  Note  that  nodes  e  U  may  be  located  on  an  edge 
of  ti  in  the  partition  7^.  Suppose  there  are  l  nodes  are  on  an  (open)edge  of  a 
triangle  ti,  we  divide  the  triangle  ti  into  l  +  1  smaller  triangles  by  connecting 
these  l  nodes  and  the  vertex  opposite  to  the  edge,  shown  in  Fig.  3.1(c).  If 
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this  edge  is  shared  by  a  pair  of  triangles  U  and  tj,  we  divide  each  of  these  two 
triangles  into  l  + 1  smaller  triangles.  Applying  this  process  to  each  edge  of  the 
triangles  will  results  in  a  triangular  partition  Tm  satisfying  (T.l)  and  (T.2). 

The  triangular  partition  Tm  of  the  cell  Qm  can  be  individually  carried  out, 
a  combination  of  the  triangular  partitions  Tm  for  all  me  Z2  forms  a  partition 
T  of  R2,  and  VT  =  UmeZ2  Km  =  UkeXfXk  where  Vj  is  the  set  of  all  vertices  of 

r. 


Based  on  such  a  partition  T,  we  can  extend  a  grid  function  u  on  Q  to  a 
function  u(x )  for  x  €  R2  by  a  linear  interpolation.  Let  4>i{x)  be  a  linear 
function  in  ti  such  that  <j>i(x^n,K^)  =  un>K  at  all  vertices  of  ti,  and  let  ipm(x)  be 
a  piecewise  linear  function  in  Qm  such  that  ipm(x)  =  (f>i(x)  in  ti,  1  <  i  <  T. 
Then,  there  holds 

|'0m(-*')  l/fi(Qm)  =  Sl<j<T  l<Ai(a')llirl(ti) 

—  u  Z^tiCQm  2-sx (.»'*) |^(n,/c)  _  x(l, A)  |2 

—  C^ticQm  E*(».«)^(W)e{j  \Un,K  —  Ui,x\2. 

Let  u(x)  —  ipm  (s)  in  Qm  for  all  m  <E  Z2.  Then,  u(x)  is  continuous  and 
piecewise  linear  function  in  R2,  and 

\un,K  ~  ul,x\2-  (3.19) 

m£Z2  tiCQm  $(n>*) ,#(*>*) Eti 

Note  that  if  the  vertices  a;(n>K)  and  x^1’^  of  ti  in  Qm  are  not  connected,  due 
to  (C'. 2)  and  Proposition  3.1,  they  are  linked  by  the  shortest  chain  Ln>Ktit\  : 
x(ni,\i)  —  x{n,n)  x(n2,x2) . . .  — >  x(na,x,)  _  x(i, x)  w^h  s  uniformly  bounded, 
where  the  node  a;K>Ai)  is  connected  to  the  node  for  1  <  j  <  s  —  1. 

Hence,  there  holds 

|wn,K  —  ul,X |  5:  C  ^  ^  \urij,Xj  ~  urij+i,\j+i  | 

1<J<S-1 

which  with  (3.19)  implies  that 

EE  E  |^m,/c  ^n+m,x\  •  (3.20) 

meZ 2  KSA/’m 


Theorem  3.8  Let  u  be  a  grid  function  on  the  lattice  Q,  and  let  u  be  the 
extension  of  u  by  a  linear  interpolation,  described  as  above.  Then  \u\h^{q)  — 
|^|ff1(R2)>  and  ||«||l2ct(^)  —  ||w|U20.(r2),  i-e.  there  are  two  positive  constants  C\ 
and  C2  independent  of  u  and  u  such  that 

C\\u\ij-L(Q)  <  |u|ffi(R2)  <  C2\u\ffl(g)  (3.21) 


Ci\\u\\Li^g)  <  \\u\\li„(R2)  <  CiWuWli^g) 


and 


(3.22) 
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Proof.  Since  u(x)  is  a  piecewise  linear  function  interpolating  the  grid  function 
u  at  each  node  x(m'K\  and  for  x  €  Qm, 

(1  +*’)■' log^l  +  |*|)  Si  (1  +  |m|2)-log^(l  +  |m|), 

it  holds  that 

ll«Wlll>,„(«„)  =  E  (1  +  Myi°«2<,(i  +  M)K.I2 

a ;(«.»)  eJTm 

which  implies  (3.22). 

The  second  inequality  of  (3.21)  follows  from  (3.20).  It  suffices  to  show  the 
first  inequality  of  (3.16).  Suppose  that  x^m,K^  and  x(n+m,x^  are  connected.  If 
x(m,K)  an(j  x(n+m,\)  are  jn  ceu  Qm,  it  is  easy  to  see  that 

|wm,K  —  ^n+m,A|  ^  C'lV’ro  Iff1  (£?,>»)•  (3.23) 


We  next  consider  two  connected  nodes  x^m,K^  and  a;(n+m>A)  which  are  located 
in  different  cells  Qm  and  Qn.  Let  Qnj,  1  <  j  <  J  be  a  sequence  of  cells  with 
Qm  =  Qm  and  Qnj  =  Qn  such  that  Qnj  is  neighboring  to  Qnj+1 .  Due  to  the 
assumption  (C".3),  J  is  uniformly  bounded.  Select  a  common  vertex  xK’-b)  0f 
the  cell  Qn .  and  Qnj+ x,  1  <j<J  —  1.  Therefore,  we  have 

^n+m,A|  ^  |  d"  I'^m+rajA  ^nj,Aj| 

(3.24) 

d"  \Unj,Xj  ~  unj+i,Xj+i  I 

Since  and  xK'+iAj+i)  are  in  the  same  cell  Qnj+1 ,  we  have  for  1  <  j  < 

J- 2 

\unj, Xj  -  Unj+1,xj+1 1  <  C\i)nj+l  \m(Qnj+1),  (3.25a) 

Similarly,  there  hold 

|^m,K  —  ^ni,Ai|  <  (3.25 b) 

and 

l^m+n,  A  ^nj_i,Aj_i  I  <  (3.25c) 

A  combination  of  (3.24)  and  (3.25)  leads  to 

|^m,K  ^n+m,X  \  —  C  7!  \lpnj\HHQnj)-  (3-26) 

1  <j<J 

The  first  inequality  of  (3.21)  follows  easily  from  (3.23)  and  (3.26).  □ 

Lemma  3.9  For  u  £  Hl(Q)  there  exists  a  constant  a  such  that 

llu  ~  a\\ L2_t  _X(6T)  <  C\u\Hi(g)  (3.27) 


Proof  Let  u(x)  be  the  extension  described  above  and  a  —  fr  u(x)dx  where 
T  =  {x  6  R2  |  |x|  =  2}.  By  Theorem  A.l  there  holds 

f  \u-a\2dx+  [  -^r-^T-dx  <  C\u\2hi{R2) 

Js  J sc  x|2log*|x  v  ’ 
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where  S  —  {a:  €  R2  |  |rc|  <  2} ,  Sc  =  R2  \  S.  This  estimation  and  Theorem  3.8 
lead  to  (3.27).  □ 

We  are  now  able  to  prove  Theorem  3.7. 

Proof  of  Theorem  3.7  Let  v(x)  be  the  extension  of  v,  and  a  —  frv(x)dx 
where  T  =  {a;  G  R2  |  |a;|  =  2}.  Since  /*  =  we  have  Theorem  3.8 
and  Lemma  3.9 

I  (Mol  =  \(f,v-a)6\ 

<  C\\f\\L2hi(g)\\v-a\\Llh_iie) 

<  C,||/llL?il(a)P  ~  _i(r2) 

<  c,II/IIl?>1(s)I^Ihi(r-2) 

<  CWfWLi^lvlmiG)- 

□ 


We  are  able  to  address  the  existence  and  uniqueness  of  the  solution  of  the 
problem  (3.1)  with  d  =  2  in  the  energy  space  E{Q). 

Theorem  3.10  If  /  €  L\  i(Q)  and  /*  =  then  the  problem  (3.1)  with 
d  =  2  has  a  solution  u  €  E(Q),  and  the 

Ml pm  <  C|l/llLi,l(e)-  (3.28) 

The  solution  is  unique  in  E{Q )  up  to  a  constant  for  s  =  1  and  a  rigid  body 
motion  for  s  —  d  =  2. 

Proof  By  the  property  of  the  bilinear  form  B 

|B(m,u)|  <  C'|MU(a)IM|£;(0) 


and 

B(u,u)  =  \\u\\2E(g), 

B  is  continuous.  Due  to  Theorem  3.7,  /  6  L\ tl(Q)  with  YlkeN  f*  =  ^  defines 
a  linear  functional  on  E(Q),  and 

E{v)  =  l(/>v)d  <  c\\f\\Llu_^G)\v\m{g)  <  C'll/llLi1>_1(a)lli;IU(e)- 

If  a  function  with  zero  strain  energy  is  regarded  as  ’’zero”  in  E(Q),  the  energy 
space  E(Q)  is  a  Hilbert  space.  By  Lax-Milgram  Theorem,  there  exists  a  unique 
solution  u  €  E(Q)  such  that 


IMIe(C0  <  CWfWq^g)- 


Hence  (3.23)  holds.  Owing  to  the  rigidity  of  the  lattices,  a  function  with 
zero  energy  is  a  constant  for  s  =  1  and  a  rigid  body  motion  for  s  =  d  =  2. 
Therefore,  the  solution  of  the  problem  (3.8)  in  two  dimensions  is  unique  up  to 
a  constant  for  s  =  1  or  a  rigid  motion  for  s  =  d  =  2.  □ 


3.3.3  Lattice  problems  without  absolute  terms  in  three  dimensions 
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For  the  existence  and  uniqueness  of  solutions  on  lattice  Q  in  R3,  we  have  to 
establish  the  extension  of  grid  functions  on  three  dimensional  lattice  Q  to  whole 
space  R3. 

As  in  two  dimensions,  let  Qm  be  the  closure  of  Qm,  and  let  Km  be  a  set  of  all 
nodes  in  Qrn.  By  K^,K^,  and  Klm,  we  denote  the  subsets  of  Km  for  nodes 
at  vertices,  on  the  edges  (not  including  the  vertices),  on  faces  (not  including 
nodes  on  the  edges  and  at  vertices),  and  in  the  interior  of  Qm,  respectively. 

By  T  =  {U,  1  <  *  <  T}  we  denote  a  tetrahedral  partition  of  Qm  satisfying 
the  conditions  : 

(T.3)  Vj-m  =  Km-,  where  Vfm  denotes  a  set  of  all  vertices  of  the  partition  T; 
(T.4)  U  fi  tj  for  i  j  is  a  vertex,  or  an  edge  of  ti,  or  a  face  of  U,  or  empty. 


(a)  Initial  partition  of  cell  Qmwith  6  simplices 


(b)  Partition  of  a  simplex  (c)  Partition  of  simplex 

t .  with  a  node  in  the  interior  t .  with  nodes  on  a  face  Fy 


Fig.  3.2  Tetrahedral  partition  Tm  of  cell  Qm 

Each  ti  is  a  simplex  with  faces  F^,  1  <  j  <  4.  The  construction  of  such  a 
partition  can  be  started  with  an  initial  partition  7^  of  Qm  such  that  (T.4)  hold 
and  Vro  =  K^,  shown  in  Fig.  3.2(a).  For  the  partition  7^  there  may  be  some 
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nodes  in  the  interior  of  simplices  Us.  If  a  node  x^n,K^  €  is  in  the  interior  of 
U,  we  divide  U  into  four  smaller  simplices  by  connecting  x^n,K^  to  four  vertices 
of  U,  shown  in  Fig.  3.2(b).  Repeating  the  process  for  each  node  in  the  interiors 
of  all  simplices,  we  have  a  partition  7^  of  Qm  for  which  (T.4)  holds  and  no 
node  a;(n>K)  €  Km  is  located  in  the  interiors  of  all  simplices.  Note  that  nodes 
x(n,K.)  ^  y  y  may  iocated  on  the  closure  of  faces  of  Us  in  the 
partition  7^.  Suppose  there  are  several  nodes  are  on  Fy  which  is  the  closure  of 
Fij.  According  to  the  triangular  partition  of  a  cell  in  two  dimensions,  described 
in  previous  subsection,  there  is  a  triangular  partition  Tptj  =  {n,  1  <  l  <  L} 
of  such  that  (T.l)  and  (T.2)  are  satisfied.  Connecting  the  vertices  of  the 
partition  Tptj  and  the  vertex  opposite  to  the  face  Fy,  we  divide  this  simplex 
U  into  several  smaller  simplices.  If  Fy  is  shared  by  a  pair  of  simplices,  the 
division  can  be  done  in  each  of  them.  Carrying  this  division  on  each  face  of 
simplices  U  and  each  simplex  in  the  partition  7^,  we  will  obtain  a  desired 
partition  Tm  satisfying  (T.3)  and  (T.4).  A  combination  of  the  partitions  Tm 
for  all  m  €  Z3  form  a  tetrahedral  partition  T  of  R3  and  Vr  =  Umez3Km, 
where  Vj  is  the  set  of  all  vertices  of  T. 


As  in  two  dimension,  based  on  such  a  tetrahedral  partition  T  of  R3,  we  can 
extend  a  grid  function  u  defined  on  three  dimensional  lattice  Q  to  a  function 
u(x)  for  x  G  R?  by  a  linear  interpolation.  Let  (f>i(x)  be  a  linear  function 
in  U  which  interpolates  u  at  the  vertices  of  U,  and  let  ipm(x)  =  </>i(x)  for 
x  €  U,  1  <  i  <  T,  which  is  a  piecewise  linear  and  continuous  function  in  Qm, 
and 


HUQm) 


<cE  E 

tiCQm  *<»•">,* V’Veti 


\'Ujn,K  Ul,\\ 

|x(n,AC)  —  |2 


<  c  k,k-ui,x\2. 

tiCQm 

Let  u(x)  =  ipm(x)  for  x  e  Qm,  m  e  Z3.  Then,  u{x)  is  a  continuous  and 
piecewise  linear  function  in  R3,  and 

<  C  £  £  £  \un,K  —  ul,\\2- 

m£Z3  tiCQm  x(n.«),o ;(!-A)eti 


Note  that  vertices  of  a  simplex  U  may  not  be  connected.  Arguing  as  in  two 
dimensions  for  (3.15),  we  have 

£1  I um,K  ~  Um+n,\\2  5;  C\u\Hi^gy 

m£Z3  neK  (n,\)€BK 

The  arguments  for  the  equivalence  between  norms  of  u  and  its  extension 
u(x )  in  two  dimensions  can  be  carried  out  in  the  three  dimensions.  Hence  we 
have  the  following  theorem  which  is  parallel  to  Theorem  3.8. 

Theorem  3.11  Let  u  be  a  grid  function  on  a  lattice  Q,  and  let  u  be  the  exten¬ 
sion  of  u  by  linear  interpolation,  described  as  above.  Then,  \u\H^(g)  —  |w|ri(r3) 
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and  \u\li  ^g)  —  \u\lijr3),  i.e.  there  are  two  positive  constants  independent  of 
u  and  u  such  that 

Ci\u\Hi(g)  <  |w|ffi(RS)  <  (3.29) 

and 

Ci\u\Ll^g)  <  |u|l2>(t(R3)  <  C2\u\Lug).  (3.30) 

Lemma  3.12  For  v  G  E(G )  there  exists  a  constant  a  such  that 

llu  _  a\\ Ll^g)  <  C\v\m(g)  (3.26) 

Proof.  Let  v  be  the  extension  of  v  described  above,  and  let 

a  =  lim  f  v(r,9,<f>)dS. 
r^°°  p|  Js 

where  S  is  the  unit  sphere  centered  at  the  origin.  By  Lemma  A.  10  the  above 
limit  exists.  Due  to  Theorem  3.11  and  Lemma  A.  12,  we  have 

lb-all  L>v)<cf  dx<cf  \Vvfdx  <  C\v\HHg). 

J  R3  r  J  R3 

□ 


Theorem  3.13  If  /  €  L\{G ),  and  Yketffk  =  °> then  for  any  v  G  E(9)i 


\(Mg\ 


Y^fkVk 

fcGJV 


<  C\\f\\Ll(g)\v\m(g) 


(3.31) 


Proof.  Let  v  be  the  extension  of  v,  and  let 

ex  —  lim  [  v(r,  6,  <j>)dS. 
r-K>o  \S\  Js 

as  in  previous  lemma.  Then  by  Cauchy  inequality  and  Lemma  3.5  and  Lemma 
3.9 

\(f,v)g\  =  \(f,v-a)g\<C\\f\\Liie)\\v-a\\Liiig) 

<  C\\f\\L2{g)\\v  -  Q'||£,2^i(K3)  <  C\\f\\L2{g)\v\Hl(R3) 

<  C\\f\\Ll(g)\v\HHg)- 

□ 

and  due  to  lemma  2.1,  its  (energy)  norm  is  equivalent  to  the  semi-norm  of 
H\Q).  Note  that  VT  —  span{ei,Ti+3, 1  <  i  <  3}  for  s  =  d  =  3,  and  Vr  =  Vc 
denotes  the  set  of  constant  functions  on  Q  for  s  =  1,  which  are  defined  in 
(2.16). 

Theorem  3.14  If  /  G  Lj(G)  and  YkeM  /*  =  °>  then  the  problem  (3.8)  with 
d  =  3  has  a  solution  u  G  E(G),  and  the 

\u\m(G)  ^  c'll/lli?,1(e)-  (3-3) 

The  solution  is  unique  in  E(G)  up  to  a  constant  for  s  =  1  and  a  rigid  body 
motion  for  s  =  d  —  3. 
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Proof.  It  is  shown  that 

\B{u,v)\  <  C'||m||e(0)||u||js(0) 

and 

B(u,u)  =  \\u\\%{g). 

By  Theorem  3.13,  f  €  L\{Q)  with  Y^kert  h  =  0  defines  a  linear  functional 
over  E(Q),  and 

\f(v)\  =  l(/>v)al  <  cil/IU^Mtf1^)  ^  c,ll/llz,f(e)llwIU(e)- 

Let  functions  with  zero  energy  be  regarded  as  a  ’’zero”  element  in  the  energy 
space  E(Q),  then  E{G )  is  a  Hilbert  space.  By  Lax-Milgram  Theorem,  there 
exists  a  solution  u  E  E(Q)  such  that 

Mffi(e)  =  IMU(s)  <  C\\f\\ L2ltl(g)-  ■ 

By  the  rigidity  of  the  lattices  and  Proposition  2.1,  a  function  with  zero  energy 
is  a  constant  grid  function  for  s  =  1  and  a  rigid  body  motion  for  s  =  d  =  3  on 
Q.  Therefore  the  solution  of  the  problem  (3.8)  is  unique  up  to  a  constant  for 
s  =  1  and  a  rigid  body  motion  for  s  =  d  =  3.  □ 

4.  UNSTRUCTURED  LATTICES  IN  HALF  SPACES 

4.1  Setting  of  unstructured  lattices  in  half  spaces  R+ 

For  a  lattice  Q  in  a  upper-half  space  R+  =  {#  =  (x',  xd)  G  Rd  |  xa  >  0},  the 
conditions  (C'.1)-(C".3)  on  connectivity  and  the  assumption  on  the  rigidity  are 
valid.  The  notations  in  previous  section  will  be  adopted  here.  In  addition  we 
need  to  precisely  characterize  unique  features  of  lattices  in  an  half  space  Let 
Td  =  {x  =  (x',xd)  G  Rd  |  xd  =  0}  be  the  hyperplane.  For  d  =  2,  T2  is  the 
real  line  and  it  is  consisting  of  rods  and  nodes.  For  d  =  3,  T2  is  the  x\  —  x-2 
plane  and  is  consisting  of  polygons  for  which  each  edge  is  a  rod  and  these 
polygons  are  faces  of  the  polyhedral  cells.  By  Kb  we  denote  the  subset  of  K 
in  which  nodes  are  on  the  hyperplane  and  K°  =  K  \  Kb.  Kb  is  referred  as 
the  boundary  of  the  lattice  Q  By  Mb  we  denote  the  subset  of  M  such  that 
Xk  G  Kb  for  k  E  Afb,  and  Af°  =  M  \  Nb. 

By  Kbm  we  denote  a  subset  of  Km  in  which  the  nodes  are  on  the  hyperplane 
Td,  and  =  Km  \  K^.  Further,  we  denote  the  subset  of  by  for  the 
indices  of  nodes  in  and  =  Afm  \  Af^.  We  may  set  md  =  0  for  those 
cells  Qm,m  E  Zd  connected  to  Td.  Therefore,  Kbm  and  Af^  are  not  empty  for 
m  E  md  =  0,  and  =  Km  and  =  Km  for  all  m  E  Z+,  md  >  0. 

4.2.  Boundary  Value  Problems  We  consider  the  boundary  value  problem 

_  V  +  C kuk  =  fk,  VA c  G  A f°  (4.1  a) 

eeKk  'X£  Xk' 

with  boundary  condition 

Uk  =  0  for  k  E  Afb  (4-16) 


(a)  A  unstructured  lattice  in  half  space  R;  (b)  An  interior  cell 

Fig.  4.1  An  unstructured  lattice  in  half  space  R/j_ 


where 

E<M>  =  AE{k’e)  —  —  ^  —  Xk^T 

\xi-xk |2  \xe-xk\2 

and 

Cfc  =  diag(Ck,  Cl  ...,Ckd),  C?>0,l<l<d. 


Let  Hl(Q)  and  L2{Q)  be  the  discrete  Sobolev  spaces  with  the  norms 

IMIW)  =  E  H2  <4-2> 

fee  A"0 


and 

IMIffi(S)  =  Mh!(£0  +  IMIi2(£?) 

where  |u|i/i(e)  is  the  semi-norm, 


m  h^g) — y, 

kerf°ieKk 


|a*  -^fcl2’ 


(4.3a) 


(4.36) 


By  Hq(Q),  we  denote  the  subspace  of  Hl{Q),  in  which  functions  vanish  at  the 
boundary  F^  =  Kb. 

The  corresponding  variational  equation 

B(u,  v )  =  F(v),  Vv  G  HlQ)  (4.4) 


are  associated  with 


(4.5a) 


is  a  bilinear  form  on  Hq(Q)  x  Hq(Q),  and 

F(v)=  E  </*.»»>• 


kerf0 


(4.56) 
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is  a  linear  functional  on  Hq{Q). 

The  energy  of  the  lattice  Q  is  defined  as 


G(u)  =  B(  u,«)=  ^{E<5e(W) 
keM°  eefck 


ll£  Uk  U£  Uk 

\xt-xky  \xe-xk\ 


)  +  (C  kuk,uk)} 


(4.6) 


As  analogue  of  the  Theorem  3.3  for  unstructured  lattice  problem,  we  have 
the  existence  and  uniqueness  for  the  boundary  value  problems  (4.1). 

Theorem  4.1  If  Ck  ^  0  in  the  sense  that  there  exists  a  node  xkm  =  x^m,Km^ 
in  each  cell  such  that  Ckm  ^  0,  the  problem  (3.1)  has  a  unique  solution  u  G 

HUG). 

4.3  Problems  Without  Absolute  Terms 

There  is  no  essential  difference  to  address  the  existence  and  uniqueness  of 
solution  for  problem  (3.6)  in  entire  spaces  and  the  boundary  value  problem 
(4.1)  in  half  spaces  if  the  problems  associated  with  absolute  terms.  We  should 
focus  the  treatment  for  problem  without  absolute  terms  in  half  space,  which 
is  essentially  different  from  those  for  the  problem  in  entire  space.  Without 
the  absolute  term  the  boundary  value  problem  is  to  find  a  u  G  Eq(G)  =  {«£ 
E(G)  |  v  =  OonTg}  such  that 

B(u,v)  =  F(v),  VveE0(G)  (4.7) 


where  b  is  a  bilinear  form  on  E0(G)  x  Eq(G), 


B(u,v)  =  S^E(M) 

keAf°(eKk 


ue  -  Uk  ve  ~  Vk 
\xe-xk\ '  \xt>-xk 


) 


(4.8a) 


and  F  is  a  linear  functional  on  Eq(G), 


F(v)  =  (M<S  =  £</*,<’*>•  (4-96) 

k€N 

For  the  problem  (4.7),  the  energy  space  E(Q)  does  not  coincide  the  Hl(Q), 
and  E(G)  is  not  embedded  in  L2(Q).  The  “energy  norm”  ||tt||s(a)  is  equivalent 
to  the  semi-norm  of  Hl{G).  Due  to  Lemma  2.1  and  the  boundary  condition 
(4.1b)  we  have  the  uniqueness  of  the  solution  in  E0(G). 

Theorem  4.2  If  u  €  E0(G)  is  a  solution  of  the  problem  (4.1)  over  the  lattice 
Q  which  satisfies  the  rigidity  assumption,  then  it  is  unique. 

Proof  If  there  is  another  solution  v  €  E0(G),  then  ||w  —  ^||e(r,|)  =  0.  Since 
the  lattice  G  is  rigid  under  the  assumption,  by  Lemma  2.1,  u  —  v  is  a  rigid 
body  motion.  Since  ( u  —  v)  |r5=  0,  hence  u  =  v.  □ 

The  theorem  tells  the  space  E0(G)  is  a  normed  space  and  the  energy  norm 
is  a  norm  for  the  space  Eo(G).  For  the  existence  of  the  solution,  we  need  to 
find  a  weighted  space  in  which  function  /  is  linear  functional  over  Eq(G).  To 
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this  end,  we  introduce  a  weighted  space  L2(G)  with  the  norm 

iMiz,2(£) =  X3  \xk>d\2)  iwfci2 

feeA/-0 

with  a  real  number  u,  where  Xk,d  is  the  d-th  coordinate  of  the  node  Xk  which  is 
the  distance  form  the  node  Xk  to  the  hyperplane  Td.  We  shall  write  L2(G)  = 
L2(Q).  Further,  we  need  to  modify  the  Hl(Q)  by  introducing  the  space  H1(G) 
with  the  norm 

IMI/rqe)  =  {iMIll^a)  “*■ 

and  the  space 

Ho(G)  =  {«  €  H\Q)  |  um,K  =  0  on  Tg} 

To  prove  E0(G)  C  EL^G),  which  is  essential  for  the  existence  of  solution,  we 
have  to  extend  the  grid  functions  on  Q  to  continuous  functions  in  R+  by  the 
linear  interpolation  in  one,  two  and  three  dimensions.  For  d  =  l,u  G  Eq(G)  is 
extended  to  u  such  that 

u{xk)  =  Uk  for  k  €  Af.  (4-10) 

For  d  =  2, 3  the  extension  of  u  based  on  the  triangular  partition  has  been 
described  in  previous  section,  can  be  carried  out  here,  and  (4.10)  holds  at 
every  node  of  G  for  all  dimension.  Arguing  as  in  proof  of  Theorem  3.8,  we 
have  the  equivalence  of  norms  of  u  and  u(x). 

Theorem  4.3  Let  u  be  a  grid  function  on  the  lattice  G,  and  let  u  be  the 
extension  of  u  by  a  linear  interpolation,  described  as  above.  Then  \u\H\^g)  = 
|u|Hi(R2),  and  IMUg^)  =  ||«||i,2)<r(R2),  i.e.  there  are  two  positive  constants  C\ 
and  C?2  independent  of  u  and  u  such  that 

C\\u\Hi(Q)  <  |u|/fi(R2)  <  C2\u\Hi(g)  (4-11) 

and 

ci\\u\\lI,AS)  -  ll«IU?,CT(R2)  <  c2\\uWl\ IAS)  (4.12) 

The  equivalence  between  norms  of  u  and  its  extension  uleads  to  a  desired 
embedding  lemma. 

Lemma  4.4  E0(G)  C  Ll1(^),and  for  u  e  E0(G),  there  holds 

MlIM)  <  C\u\m(g)  (4-12) 

with  C  independent  of  u. 

Proof  Let  u(x )  be  the  extension  of  u  by  a  linear  interpolation  described  above. 
Due  to  Theorem  4.3,  u(x)  G  Eq(R+).  By  Lemma  B.l, 

ll^llili(R+)  - 

which  together  with  Theorem  4.3  imply  (4.12)  immediately.  □ 

The  following  lemma  indicates  that  the  energy  norm  ||w||b(^)  is  a  norm  of 
the  energy  space  Eq(G)- 
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Lemma  4.5  The  space  E0(Q)  is  equivalent  to  the  space  Hq(G),  and  for  u  G 
E0(Q)  there  holds 

IMU(0  -  II«IIhi {g)  -  Hhhg)-  (4T3) 

□ 

Lemma  4.6  If  /  G  L\(Q),  then  for  any  v  G  Eq(G)  there  holds 

l</,«->l<C||/||iW)|HUW).  (4.14) 

Proof  By  Schwarz  inequality,  there  holds 

|(/,w)|  < 

Due  to  Lemma  4.4-4. 5  we  have 

K/j  v) I  <  C\\f\\Li(g)\v\H'(g)  <  C'II/IIl?(0)IMU(0)- 

□ 

The  above  lemmas  lead  to  the  existence  and  uniqueness  of  solution  for  the 
problem  (4.7). 

Theorem  4.7  For  /  G  L\{Q),  the  problem  (4.7)  for  d  =  1,2,3  has  a  unique 
solution  u  G  E0(G),  and 

IMUw)  <  Cll/llqw)-  (4-15) 

Proof  Due  to  Lemma  4.5,  we  have  for  it,  v  G  Eo(G) 

\B(u,v)\  < 

and 

B(u,u)  =  >  -^IMI#i(g)* 

By  Lemma  4.6,  /  defines  a  linear  functional  F(v )  on  E0(Q),  and  for  v  G  E0(Q) 
l-F(u)l  =  \(fiv) I  <  CWfWLi^lMEig)  <  C'll/llLi1(e)llvllffi(6)- 
By  Lax-Milgram  Theorem,  there  exists  a  unique  solution  u  G  Hq(Q)  such  that 

IMlirqc;)  ^  ^ll/IUlde) 

which  with  Lemma  4.6  leads  to  the  assertion  of  the  theorem.  □ 

Remark  4-1  The  boundary  T#  can  be  a  very  general  curve,  and  the  function 
xd  =  4>(x')  is  required  to  be  a  piecewise  continuous  function  and  allow  to 
approach  oo  as  |a/|  — >  oo,  e.g.  £2  =  x2  for  d  =  2.  The  weight  (1  +  xdy  may 
be  modified  to  (1  +  d(x)2)u  where  d(x)  is  the  distance  from  x  to  the  boundary 
r*,  for  the  detalis  of  such  a  modification  we  refer  to  [2]. 

CONCLUSION 

The  lemmas  and  theorems  in  previous  sections  indicates  that  a  mathematical 
framework  for  problems  of  unbounded  lattices  has  been  established.  In  this 
framework  we  are  able  to  prove  the  existence  and  uniqueness  of  solution  for  the 
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problems  of  unstructured  lattice  in  entire  and  half  spaces.  This  is  a  significant 
progress  in  research  of  lattice  problems. 

We  have  derived  the  appropriate  function  spaces  which  the  data  /  (e.g. 
external  force)  belongs  to,  and  proved  the  existence  of  solutions  for  the  lattice 
problems  without  absolute  terms  in  entire  and  half  spaces  if  the  date  /  is 
given  in  these  function  spaces.  The  function  spaces  and  approaches  to  prove 
the  existence  theorems  are  quite  different  in  one, two  and  three  dimension. 
Meanwhile,  because  of  the  Dirichlet  boundary  condition  for  lattice  problems 
without  absolute  terms  in  half  spaces,  the  function  spaces  and  approaches  to 
prove  the  existence  theorems  are  the  same  or  very  similar  in  one,  two  and 
three  dimension.  This  reflects  the  fundamental  differences  between  problems 
without  boundary  and  the  boundary  value  problems.  The  results  reported 
here  and  in  forth  coming  papers  could  be  further  extended  and  generalized 
as  indicated  in  various  remarks  and  comments  in  previous  sections.  Although 
the  theorems  are  proved  for  truss  problems,  it  can  be  utilized  or  generalized  to 
problems  of  general  lattices  such  as  plates  shells,  and  three  dimensional  solid. 

Based  on  the  progresses  we  have  made  in  the  past  years,  our  research  on 
lattices  should  be  further  carried  on.  The  focus  of  the  research  at  next  stage 
will  be  the  analysis  of  complicated  models  and  design  of  effective  computation. 
A  direct  extension  of  our  research  results  is  the  investigation  on  the  existence 
and  uniqueness  of  solutions  for  general  lattice  problems  without  absolute  terms 
in  unbounded  domains,  associate  with  various  structures  such  as  plate,  shell 
and  3-dimensional  solid.  This  is  an  important  issue  and  has  been  an  not-  well- 
unanswered  problems  for  last  2  decades.  The  extension  of  grid  functions  by 
a  linear  interpolation  for  rigid  and  non-triangular  or  non-tetrahedral  lattices 
implies  an  effective  numerical  approach  for  unstructured  lattices,  i.e.  multigrid 
method.  To  this  end  we  shall  define  partial  differential  equations (PDE)  asso¬ 
ciate  with  boundary  conditions  or  no  boundary  conditions,  which  is  equivalent 
to  corresponding  lattice  problems.  This  equivalent  PDE  is  not  homogenized 
one  because  homogenization  is  not  applicable  to  the  unstructured  lattice,  and 
can  be  solved  on  fine  and  coarse  grids.  The  discretization  of  the  PDE  on  the 
grid  of  lattices  is  the  original  problem  of  lattice.  The  approach  has  been  inves¬ 
tigated  for  unstructured  lattices  in  bounded  domains  [4],  it  could  be  carried 
for  unbounded  lattices. 

If  the  scale  of  cells  of  lattices  is  very  small,  we  have  to  deal  with  a  multi-scale 
problem  which  is  one  of  the  most  popular  topics  of  modeling  and  computing  in 
modern  applied  and  computational  mathematics  and  engineering.  The  frame¬ 
work  we  have  built  for  the  problem  without  absolute  terms  will  help  us  to 
analyze  various  aspects  of  such  multi-scale  problems. 
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APPENDIX  A 

In  this  appendix  we  will  address  the  existence  and  uniqueness  of  solution  of 
the  Dirichlet  problem  of  Laplace  equation  in  Rd,  d  =  2, 3.  We  are  interested  in 
this  problem  because  it  is  similar  to  the  problem  (3.1)  of  unstructured  lattices 
in  entire  spaces  without  absolute  term.  In  particular,  some  embedding  results 
play  a  key  role  for  proving  the  existence  and  uniqueness  of  solution  for  such 
lattice  problems.  We  seek  a  solution  u  e  H(Rd)  of  the  variational  problem  : 

B(u,  v )  =  F{v)  Vv  G  H{Rd)  (Ala) 


(Ale) 


B(u,  v )  =  F(v)  Vv  €  H{Rd)  (Ala) 

with 

B(u,v)  =  [  Vu-Vvdx,  u,v  e  H(Rd)  (A lb) 

J  Rd 

and 

F(v)  =  f  fvdx,  v  e  H(Rd)  (Ale) 

J  Rd 

where  d  =  2,3,  and 

H(Rd)  =  {u\ue  Hlocal(Rd),  MHi(Rd)  <  oo} 

Here  |u|i/i(Rd)  is  the  semi-norm  of  the  Sobolev  space  H1^)  involving  only 
the  first  derivatives. 

A.l  In  two  dimensions 

Let  H1  (R2)  be  a  closure  of  C00  functions  with  the  norm 


IMHw  -  Js  I“I  -  Jsc  r2 log2  '  JK 2 1  v  v^'“' 

where  (r,  $)  are  the  polar  coordinates,  the  disc  S  =  {x  G  R2|r  =  \x\  <  2}  and 
Sc  —  K2\  S.  By  L2  (R2)  we  denote  a  weighted  space  with  the  norm 


\u\2dx  + 


[  |V«|a 

J  R2 


(A2) 


IMII2i0(r2)=/  M 2dx+  [  |u|2r2l/  log2*7  rdx.  (A3) 

J  S  J  Sc 


Then  we  have 


We  further  introduce 


H0(R2)  =  |tt  €  H{ R2)  |  J  uds  =  0, 
Hr(Sc)  =  jw  |  J  \Vu\2dx  <  a,u\r= 


-} 


where  V  —  dS,  and  we  define  a  quotient  space 


with  the  norm 


H\ R2)  =  H\R2)/Pq 


|«||  *  =  inf  || u  -  a||jyi(R2) 

hhr2)  aepo 


(A. 4a) 


(>1.46) 
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where  P0  is  the  set  of  all  real  numbers.  Since  P0  C  H1( R2),  the  quotient  space 
is  well  defined.  For  the  existence  and  uniqueness  of  the  solution  of  the  problem 
(A.l)  with  d  =  2,  we  need  to  establish  the  equivalence  between  ||w||  ~  and 

|u|hi(R2)  and  to  show  that  F(v)  in  (Ale)  defines  a  linear  functional  over 
Hl(  R2). 

Theorem  A.l  If  u  €  Ho  (R2),  then 

[  \u\2dx  +  [  -^^-dxKC  [  \Vu\2dx  (A5) 

Js  JSc  r2  log  r  Jr 2 

with  constant  C  independent  of  u. 

Theorem  A.l  is  a  major  theorem  in  this  section,  we  need  several  lemmas  to 
prove  it. 

Lemma  A.2  Let  u(t)  be  a  function  on  (0, 00)  satisfying  «(2)  =  0  and 


l 


°°'du 
dt 


tdt  <  00. 


Then 


J2  t  log2 1  h 


tlog2t 

Proof.  Due  to  Theorem  1.14  of  [29] 


poo  pc 

J  \u\2wdt  J 


00  du 
dt 

du 
dt 


vdt 


vdt 


where  w  =  j^~t,  v  —  t  and  CL  =  sup2<t<oo  Fi(t),  with 


Note  that 


and 


W)  =  (/ 

/OO 

wdt  - 


wdt 


s: 


dt 


1 2  t log2 1  log t 

f  v~ldt=  [  jdt  =  log<  -  log  2. 
J  2  J  2  * 


Then 


FL(t)  = 


flogt  —  log2 


log* 


<  1  for  2  <  *  <  00, 


and  (A. 6)  follows  immediately. 
Lemma  A. 3  If  u  €  Hr  (Sc),  then 


f  AAL-dX<c[  \vuf 

JSc  r2  log  r  Js c 


dx 


(A- 6) 


(A.  7) 
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Proof.  Note  that  u(r,  6)  |r-2=  0,  which  implies  by  Lemma  A. 2  that 


rdr 


with  constant  C  independent  of  8.  Integrating  with  respect  to  9  from  0  to 
we  have  (A.6). 


Lemma  A.4  Let  u  e  if1  (5),  where  H 1(S)  is  a  usual  Sobolev  space  on  a 
bounded  domain  S.  Then  the  following  norm 


IMIIhhs)  - 


is  equivalent  to  the  norm,  ||it||/fi(s). 

Proof.  Obviously  |||«|||ffi(5)  is  a  norm  to  #*(5).  Note  that 

II 2 

J  uds\  <  C  \u\2ds  <  C'||«ll^(r)  <  C'll«lli?1(‘S) 

which  implies 

||M||ffi(S)  <  C'i||u||iri(S). 

We  need  to  show  that 

IMIffi(S)  <  ^lllwlllirqs). 

If  it  is  false,  there  exists  a  sequence  Uj  <G  j  =  1,2,...  such  that 

llwjllfrqs)  =  1>  and 


\u 


7 II  lff1(S)  —  lujlff2(S)  +  J 


uds 


-4  0  as  j  — >  oo 


Since  H1  (S)  CC  L2(S),  there  exists  a  subsequence  denoted  by  %  again,  which 
is  a  Cauchy  sequence  in  L2(S).  Since  \uj\H2(s)  — >  0  as  j  — f  oo,  {uj}(^_1  is  a 
Cauchy  sequence  in  if1  (5)  as  well.  Hence  limj_»oo  Uj  =  uq  in  Hl(S).  This 
implies  that  Dau0  =  lira,-^  DaUj  =  0  for  |a|  =  1.  Therefore  u0  is  a  constant 
in  S.  Note  that 


I  <c£  («,  -  «o)2  ds  <  C\\v.j  -  «ofel(S)  o  as  j  -4  oo, 

which  leads  to 

/  Unds  =  lim  /  Ujds  =  0. 

Jr  i-+°°  J 

Hence  uq  =  0  in  S.  It  contradicts  the  fact  that  ||ito||tfi(s)  =  lim^oo 
||itj||ffi(5)  =  1.  Thus  the  lemma  is  proved. 


Lemma  A.5  If  u  €  H0  (R2) ,  u  has  the  Fourier  series  on  T: 


OO 

u(  2, 0)  =  ^2ak  cos  +  bk  sin  kd 

k= 1 


47 


y,  (°fc + k  <  c\u\2Hi(S)  <  c'i«i^i(R2) 


Proof.  Since  fr  uds  =  0,  the  coefficient  ao  =  0,  and 

OO 

U (2, 9)  =  y:  (afe  cos  kd  +  sin  A;0) . 
Then  by  the  trace  theorem  we  have 

(X) 

y~l  (afe  +  &*)  A:  «  Ilullffi/2(r)  <  Cf|lMll/ri(s) 


by  Lemma  A.4 


We  now  prove  Theorem  A.l. 

Proof  of  Theorem  A.l.  For  u  e  H0  (R2)  we  can  find  a  harmonic  function 
ui  such  that  ( u  —  «i)  |r  =  0  and  fsc  |Vwi|2  dx  <  oo.  Let  w2  =  u  —  U\.  Then 
u2  0  Hq  (Sc),  and  by  Lemma  A. 3 


u2\2  *  -  dx  <C  f  |Vu2|2  dx 

rl  log  r  Jsc 


I  0  1  2  —  |  V  | 

ySe  r2  log  r 

Since  ui(2,6)  =  u(2,9),ui  has  a  Fourier  series  on  T 

OO 

Mi  (2, 9)  —  ( ak  cos  k9  +  bk  sin  k9 ) 

fc=i 

with  ao  =  0.  Because  U\  is  harmonic  in  Sc, 


ui  (r,  9)  =  ^  (a*  cos  k9  +  bk  sin  A :9) 


(A- 9) 


L  '“l|2  ^~rrdrM  5  CE«  +  «4 


00  rdr 

^.2fc+2  T 


<  C'Zirt  +  bl)* 


by  Lemma  A. 5 

<  C'|u|^ri(S) 

Note  that  w2  |r=  0  and  A«i  =  0  in  S'0,  which  implies  that 

f  Vu\Vu2dx  =  f  u2^Ads  —  f  u2Auidx  =  0. 


(A.  10) 


Hence 


Mlr^S®)  —  +  IM2|h1(5c 
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which  together  with  (A.9)  -  (A.  10)  leads  to 

<  C  +  |w2|ffi(5c)j  (^-11) 

—  C'|m|^1(R2). 

We  next  shall  show  that 

IM|i,2(5)  <  C\u\Hi(S)  (A- 12) 

Let  Vi  be  harmonic  in  S  and  ui|r  =  u\r,  and  let  V2  =  u  —  V\.  Since  t^lr  =  0, 
by  Lemma  A.4 

IM|l2(S)  <  IM|i?i(S)  ^  &  Mtfi(S)  •  (A.13) 

Since  V\  is  harmonic  and  fr  v^ds  —  fr  uds  =  0 

00  (r\k 

vi  (r,  0)  =  (ak  cos  kO  +  bk  sin  k6) 

fc=i 

and 

DO 

lwill2(S)  ^  C  (°*  A  h  —  C\u\ ^1(5)-  (^-14) 

k= 1 

Since  V\  is  harmonic  and  v%  vanishes  on  T,  we  have,  by  the  argument  above 
for  u\  and  «2>  that 

which  together  with  (A.13)  -  (A. 14)  leads  to  (A.12)  immediately.  A  combina¬ 
tion  of  (A.ll)  and  (A.12)  yields  (A. 5). 

Corollary  A.6  The  norm  ||u||#i(R2)  is  equivalent  to  |u|#i(R2),  and  the  spaces 
H1( R2)  and  H( R2)  are  equivalent. 

Remark  A.l.  The  weight  function  w(x )  =  1  in  S,  and  w(x)  =  r_2log_2r  in 
Sc.  Sc  excludes  the  origin  and  unit  circle.  We  may  select  others  weight,  e.g., 
w(x)  =  (1  +  r2)-1  log-2  (2  +  r)  for  all  x  G  R2.  It  is  essential  for  the  selection 
of  the  weight  that 

\w\  =  O  (|x|-2log-2  |a:|)  for  large  |rr|. 

Also  S  can  be  selected  to  any  bounded  domain  with  Lipschitz  boundary,  and 
it  is  not  necessary  to  be  the  disk  centered  at  the  origin  and  with  radius  2. 


Theorem  A.7  If  /  e  L\  x  (R2),  and  fR2  fdx  =  0,  then  for  any  v  €  H  (R2) 


/ 

Jr2 


fvdx 


<  C||/Ilif  ,(R2)lt;|£r1(R2) 


with  constant  C  independent  of  u. 


(A15) 
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*  *. 


Proof.  Let  v  £  H  (R2),  and  let  a  =  fr  vds.  Then  (v  —  a)  £  H°  (R2),  and 


/  fvdx 

= 

/  f(v  —  a)dx 

\J  R2 

7r2 

<  1 

ll/|li2)1(R2)lb-a 

by  Theorem  A.l 


<  C^II/IUf  1(r*)|v|hi(r2)- 


Theorem  A.8  If  /  G  L2  x  (R2)  and  Jr2  fdx  =  0,  then  the  problem  (A.l)  with 
d  =  2  has  a  solution  u  €  H  (R2),  and 

Mh^r2)  <  C\\f\\L2i{K2).  (A.16) 

The  solution  is  unique  up  to  a  constant. 

Proof.  Due  to  Corollary  A.6  we  have 


and 


\B(U,V)\  <  |u|Rl(R2)|w|Hl(R2)  <  C\\u\\  ~  ,  IMIffl(R) 

"  (R  ) 


\B(u,u)  \  =  |«|R1(R2)  =  |M|2~  . 

V  ifi(R2) 


By  Lemma  A.2  there  holds 

1^)1  <  C'H/||ifil(R2)l?;|ff1(R2) 

=  c||/||Lf,lP»)|Mi 


R!(R2) 


By  Lax-Milgram  lemma  there  exist  a  unique  solution  u  G  H1  (R2)  such  that 


Mi,,*,,  s  oifhum 


which  together  with  Corollary  A.6  leads  to  the  assertion  of  the  theorem. 

Corollary  A.9.  If  /  €  L20  (R2)  with  v  >  1  the  problem  (A. 15)  has  a  unique 
solution  in  H  (R2)  up  to  a  constant,  and 

IMIh(r2)  ^  c||/IIl20(r2)* 


A.2  In  three  dimensions 

In  three  dimension  we  shall  deal  with  the  existence  and  uniqueness  of  the 
solution  of  variational  problem  (A.l)  in  the  way  different  from  the  approach 
in  two  dimensions,  e.g.  we  will  not  introduce  the  quotient  spaces. 

Let 

u{r)  =  |^j  jf  u(r,  9,  <f>)dS  =  j|y  jf  u{r,  9,  <j>)sin9  d9  d<f>  (A.  17) 

where  S  denotes  the  unit  sphere,  and  (r,  9, 4>)  are  the  spheric  coordinates. 
Lemma  A. 10  If  u  €  iI(R3),  then  lim^oo  u(r)  =  A  exists. 
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Proof.  For  r  >  1, 


a(r)  =  W\Ui 


du(t,  9, 4>) 
dt 


dtdS  +  u(l) 


Let  rj,j  =  1,2,...  be  an  arbitrary  sequence  with  Hindoo tj  =  oo.  For  r,  >  r, 
we  have 

rrj  du(t,  9 ,  (j) ) 


h)-a(ri)i  =  w\LL' 
-  c(l£ 


dt 


-dtdS 


du(t,  9, 4)  |2  2 


dt 


fdtdS 


T(L 


N  i/2 

t~2dt J 

n  / 


(1  1\1/2 

-  c  \7t "  7,) 


This  implies  that  {w(rj)}?i x  is  a  Cauchy  sequence  and  that  u(rj)  converges 
to  the  same  limit  A  for  all  sequence  {r^}^  with  lim^oo  rj  =  oo.  Therefore, 
linv-^oo  u(r)  exists,  and  lim,.-**,  u(r )  =  A. 

Lemma  A. 11  Let  u(r)  be  given  in  (A. 17)  and  A  =  limr^.oo  u(r).  Then 

POO  POO 

/  \u(r)  -  A\2dr  <  C  |u'(r)|2r2dr  <  C'|m|j?i(r3)  (A.18) 

Jo  Jo 

Proof  Let  w(r )  =  u(r)  -  A.  Then  limr_>00u;(r)  =  0.  By  Hardy  inequality  330 
[19],  we  have 

POO  POO 

/  \w(r)\2dr<C  /  \w'(r)\2r2dr 

Jo  Jo 

which  is  the  first  inequality  of  (A.18),  The  second  one  follows  from 


POO  POO  P 

/  |u'(r)|2r2dr  <  / 

Jo  Jo  Js 


r*drdS. 

or 


□ 


Theorem  A.  12  If  |u|iji(R3)  <  oo,  then  there  exist  a  constant  a  such  that 

f  Z^-dx  <C  [  \Vu\2dx. 

Jr 3  r2  JR s 

Proof.Let  a  =  A  =  lim^oo  u(r).  Then 

L  s  0  (/..  T4"  *  /.,  * 


-dx  < 


r»2  1  /_o  0*2 

For  the  first  term  on  the  right  hand  side  of  (A.  19),  we  have  by  Lemma  A.  11 
|u(r)  —  a12 
lR3  r2 

For  the  second  term  we  write 

u(r,  9, 4)  -  u(r)  =  t4t  /*  ( u(r ,  9,  <j>)  -  u(r,  9',  <j>'))  sin9'  d9'  dfi 
\d\  Js 


19) 


L 


[  f  |u(r)  -  a\2drdS  <  (A.20) 

Js  Jo 


4 


t  'fc. 
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and 


u(r,  9,  <f>)  -  u(r,  9',  <f>')  =  u(r,  9,  <f>)  -  u(r ,  6',  </>)  +  u(r,  O',  4>)  -  u(r,  9', 
Note  that 

\u{r,6,(j)) -u(r,9',<t>)\2  =  |  [ 

Je> 


which  implies 


and 


p  7T 

<  C  | 

JO 


dujryr,  2 
e<  dr 
*  .du(r,  9,  <j>).2 


d9 

Mr,  9,  <f)  |2 
39 


d9 


dS 


{A.  21) 


(>4.22) 

□ 


[  Mr, M) -«(r, v_M_dx<cf  | Vu?dx. 

Jr  a  r2  ./R3 

Similarly,  it  can  be  shown  that 

[  [  ivui p*. 

/r3  ./r3 

A  combination  of  (A.19)-(A.22)  leads  to  (A.18). 

Theorem  A.  13  If  /  G  Lf(R3),  and  /r3  fdx  =  0,  then  for  any  v  €  if^R3), 

I  [  fvdx\  <  C'||/||i2(R3)|t;||ri(R3).  (>4.23) 

./R3 

Here  and  thereafter  Dr  denote  a  ball  centered  at  the  origin  with  radius  R,  and 


Proof.  Since  fR3  fdx  =  0 


/  fvdx 

J  R3 


[  f(v  —  A)dx 

J  R3 


with  A  =  limr_).oo  tttt  fs  v{r,  9,  (j))dS.  By  Lemma  A.ll, 

Pi 

I  Jr3  fvdx I  ^  C  {/r3  I/I2(l  +  r2)dx}  |/R3  ^  _(_  r2^x| 
<  C'll/llLf(R3)lv|H1(R3)- 


1/2 


□ 


The  next  theorem  addresses  the  existence  and  uniqueness  of  the  solution  of 
the  variational  equation  (A.l)  with  d  =  3  in  the  case  that  fR3  fdx  =  0. 

Theorem  A.14  If  /  €  L2(R3),  and  fRS  fdx  —  0,  then  the  variational  problem 
(A.l)  with  d  =  3  has  a  solution  u  €  iL(R3), 

M//i(R3)  —  ^II/IIl?(R3)i  (A.  24) 

and  the  solution  is  unique  up  to  a  constant. 


* 
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Proof.  It  is  easy  to  see  that  for  u,v  G  H (R3) 

\B{u,v)\  <  C|«|jsri(R8)|v|iri(R8), 

\B(u,u)\  =  |«||ri(R») 

and  by  Theorem  A.ll,  /  is  a  linear  functional  on  H( R3),  and 

|  [  fvdx |  <  C'||/||£,2(R3)|,y||ri(R3). 

JR3 

If  functions  v  in  H( R3)  with  M^1(r3)  =  0  is  regarded  as  a  ’’zero”  element, 
then  the  space  tf^R3)  is  a  Hilbert  space  with  ||u||h(rs)  =  |«|ri(r8).  By  Lax- 
Milgram  Theorem,  there  exist  a  unique  solution  u  G  H( R3),  and 

IMIffCR3)  <  C'||/|Il?(r3)- 

which  leads  to  the  assertion  of  the  theorem.  □ 


In  three  dimensions  we  noticed  that  there  are  functions  «o  €  tf(R3)  satisfies 
the  equation 

-Au0  =  /o- 

with  /R 3  f0dx  #  0,  and  «0(«)  =  O(i)  as  \x\  ->•  oo.  For  example, 


Uo  =  \ 


and 


CM 

CO 

in  D\ 

I  107T 

3_(3 
<  8tt  2 

i— 1 1 

1 

in  D{ 

'  3 

in  D\ 

A«o  =  i 

47 r 

>  o 

in  D{ 

where  D\  is  a  unit  ball  and  D\  —  R3  \  D\.  Note  that  Jr3  f0dx  =  1  and  /0  has 
compact  support  D\,  and  luoli/qR3)  <  oo.  Hence  in  the  case  that  fRS  fdx  ^  0 
we  seek  a  weak  solution  u  —  coito  +  w  G  H( R3)  such  that  w  satisfies  the 
variational  equation  : 

B(w,  v)  =  f  (f-  c0fo)vdx  Vu  G  H( R3)  (-4.25) 

Jr3 

with  co  =  /R 3  fdx.  Since  z  =  (/  —  c0/o)  G  Lf(R3)  and  Jr3  zdx  =  0,  by 
Theorem  A.  13,  w  exists,  and 

IMIff(R3)  ^  C\\z\\L2{r3)  <  C(\c0\  +  ||/||if(R3))- 


Therefore  we  have  a  corollary  on  the  existence  and  uniqueness  of  the  solution 
for  variational  equation  (A.l)  in  three  dimensions. 


1  A 
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Corollary  A. 15  If  /  €  Z^(R3),  and  Jr3  fdx  =  Co  7^  0,  then  there  exists  a 
solution  u  €  tf(R3)  such  that  u  =  c0u0  +  w,  where  w  €  (R3)  is  the  solution 

of  the  variational  equation  (A.  25), 

Mffl(R3)  <  C'(||/||z,2(R3)  +  |  f  fdx\). 

J  R3 

and  the  solution  is  unique  up  to  a  constant. 


4 
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APPENDIX  B 


In  this  appendix  we  will  address  the  existence  and  uniqueness  of  solution 
of  the  Dirichlet  problem  of  Laplace  equation  in  half  space  or  in  domains  of 
half-space  type.  We  are  interested  in  this  problem  because  it  is  similar  to 
the  problem  (3.1)  of  unbounded  lattices  with  boundary  condition  and  without 
absolute  term.  In  particular,  an  embedding  result  contained  in  Lemma  A.l 
play  a  key  role  for  proving  the  existence  and  uniqueness  of  solution  for  such 
lattice  problems. 

Let  5+  —  {z  =  (x',xd)  G  Rd  |  xd  >  0}  be  the  upper  half  space,  and  let 
Td  =  {x  =  (x',xd)  G  Rd  |  xd  =  0}  be  its  boundary.  By  5X(5^)  we  denote  the 
Sobolev  space  over  5+,  and  we  introduce  the  energy  space  : 

5(5+)  =  ju  |  u  G  Hlocal(R+),  |w|jyi(#d)  <  ooj 

where  \u\Hi^Rd^  is  the  semi- norm  of  the  Sobolev  space  i5(5+)  involving  only 
the  first  derivatives.  5(5+)  and  5X(5+)  are  not  equivalent,  and  the  energy 
norm  is  equal  to  the  semi-norm  \u\HyRd_y  By  5o(5+)  we  denote  the 

subspace  of  5(5+)  in  which  functions  vanish  on  Td. 

We  now  address  the  existence  and  uniqueness  of  the  solution  of  the  varia¬ 
tional  problem  :  seek  u  G  5o(5+),  d  =  2,3  such  that 

B(u,v)  =  F(v),  VveE0(Rd+)  (5.1a) 


with 


and 


(5.1c) 


B(u,v)=  f  Vu-Vvdx,  u,v  £  E0(R+)  (5.16) 

Jr% 

F(v)  —  f  fvdx,  v  G  E0(Rd ) 

Jr l 

where  E0{Rd+)  =  {«G  E(Rd+ )  \  u  =  0  on  rrf},  Td  =  {x  =  (x',xd)  G  Rd  \  xd  = 
0}  is  the  hyperplane,  and 

E(Rd+)  =  {u\ue  H}ocal(Rd+),  <  oo} 

Here  \u\Hi^R^  is  the  semi- norm  of  the  Sobolev  space  Hl(R+)  involving  only  the 
first  derivatives,  and  R+ =  {x  =  (x1,  xd)  G  Rd  \  xd  >  0}  is  the  upper-half  space. 
5(5+)  is  called  the  energy  space  with  energy  norm  |m|e(e^.)  =  5(u,  v)1!2  which 
is  equivalent  to  the  semi-norm  \u\HyRd_y  5o(5+)  is  its  sub  space  of  functions 
vanishing  at  the  boundary  which  is  not  equivalent  to  the  Sobolev  51(5+), 
and  is  not  embedded  in  L2(5+).  Therefore,  the  problem  (3.1)  may  have  no 
solution  in  50(5+),  e.g.  if  /  G  L2(5+).  Hence  we  have  to  find  a  space  for 
/  which  has  stronger  topology  than  L2(5+)  for  existence  of  solutions  for  the 
problem  (B.l).  To  this  end,  we  introduce  a  weighted  space  L2(5+)  furnished 
with  a  norm 

(5.2) 


INI \  M2(1  + 
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We  further  modify  the  Hl(R$.)  by  introducing  the  space  Hl(R%)  with  its  norm 

=  i(/Jd )  +  Rd)-  (•B-3) 

and  the  space  Hq(R^.)  =  {H1(R‘l)  \  u  =  0  on  1^}. 

We  are  interested  in  such  a  problem  because  it  is  parallel  to  the  problem  (4.1) 
of  unbounded  lattices  without  absolute  term,  and  a  embedding  result  contained 
in  Lemma  B.l  play  a  key  role  for  proving  the  existence  and  uniqueness  of 
solution  for  such  lattice  problems. 


Lemma  B.l  Eq(R%)  C  L2_1(R<1),  and  for  u  G  Eq(R%) 

with  constant  C  independent  of  u. 

Proof  For  u  G  Eq(R+),  there  holds 


u(x)  =  u(x',xd)  = 

Jo 


Xd  du(x',  t) 
dt 

Selecting  a  a  G  (0, 1),  we  have  by  Schwarz  inequality 

k*)is 


dt 


dt 


dt 


rdt 


=cxrn 


which  leads  to 


(BA) 


(B.5) 


ini*.,*,  <c/4^ri^p«)^ 
= cf*-.{r  (rh?/od 


(B.6) 


Note  that 


+ 

idu(^1 


£ 


1  — (T 


;dXd)dt 


/0  1  & 

Combining  (B.5)-(B.7)  we  have 


1  1  *  1  +  ar^' 

=  c/( 


i) , 


oo  j  5  v)  12^ 


(B.7) 


-  C\u\H\(Rd+-)  = 

□ 

As  a  consequence  of  Lemma  B.l,  we  have  the  following  corollary  which 
indicates  that  the  space  E0(R+)  is  a  normed  space  and  the  energy  norm  is  a 
norm  of  the  space  Eq(R+),  which  is  equivalent  to  the  norm  of  H1(R<^). 
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Corollary  B.2  The  space  Eq(Q)  is  equivalent  to  the  space  Hl(Rd+),  and 
||u||e(£)  is  a  norm  of  the  space  E0(Rdr).  For  u  G  Eo(R+),  there  holds 

IIuIIe(r^)  —  M#1^)  —  IMI/r1  (#{.)•  C®-^) 


Lemma  B.3  If  /  G  Z/?(i?+),  then  for  any  v  G  Eo(R+),  there  holds 

I </>«)!  <  C\\f\\q{Rd+)\\v\\2E{Rd+y  (B. 9) 

Proof  By  Schwarz  inequality,  there  holds 

\(f,v)\  <  II/L^IMIlI^) 

which  together  with  (B.4)  and  (B.8)  leads  to  (B.9)  immediately.  □ 

We  now  are  able  to  address  the  existence  and  uniqueness  of  solution  for  the 
problem  (B.l)  in  the  framework  of  the  space  Hq(R+). 

Theorem  B.4  For  /  G  L\{Rd+ ),  the  problem  (B.l)  has  a  unique  solution 
u  G  E0{R+),  and 

IMI mi)  <  c\\f\\LW.  (b. io) 

Proof  By  Corollary  B.2,  there  hold  for  u,v  G  Eo(R+) 

\B(u,v)\  < 

and 

B(U,U )  =  >  •^IImIIh1(b^)- 

Due  to  Lemma  B.3,  /  G  L\(Rd+)  defines  a  linear  functional  on  Eo(R{),  and 

|T’(u)|  <  C'||/||£,2(Hd)||u||£;(Bci_)  <  C,||/||£,2(Jid)||v||^l(Jjd). 

By  Lax-Milgram  Theorem,  the  variational  problem  (B.l)  has  a  unique  solution 
u  G  Hl{R%),  and 

—  ^ll/IU?(a)- 

Corollary  B.2  implies  the  existence  and  uniqueness  of  the  solution  in  Eq(R+) 
for  the  problem  (B.l),  and  (B.10)  follows  from  (B.8). 
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